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ABSTRACT 


The  problems  of  diffraction  of  scalar  waves  by  a  circular  aperture 
in  a  perfectly  soft  and  in  a  perfectly  rigid  infinite  planar  screen  are 
studied.   Excitation  of  the  form  cos  m{^-^^^'^' )u^^' {r,z) ,    m  =  0,1^2,  ...^ 
originating  in  the  region  z  <  0  is  assumed.   Here,  with  the  z-axis  taken 
as  the  axis  of  symmetry,  r,  ^  and  z  denote  th^e  cylindrical  coordinates 
of  an  arbitrary  space  point  and  ^   ^^'    is  an  arbitrary  phase  angle.   On 
the  shadow  side  of  the  screen,  the  solution  V   „    is  represented  as 
follows : 

V^^2=  coM^-Zj'hr^^  ^"^    J   rexprikR-'(r,z;t)]V(r,z;t)lf^^2^'")(t)dt. 

The  subscripts  1  and  2  in  this  expression  refer  to  the  solutions  for  the 
soft  and  rigid  screen  problems  respectively;  k  is  the  wave  number 


+  /2  2 

^  (r,z;t)  =  /y/r  +  (z  +  iat)  ; 


z  =  0  defines  the  plane  of  the  screen;  'a'  is  the  radius  of  the  aperture, 

(m),        (m), 
centered  at  the  origin.   The  unknown  functions  f"  [tj    and  r  itj  are  re- 
quired, respectively,  to  be  odd  and  even  analytic  functions  of  t,  to  be 
regular  in  a  neighborhood  of  the  unit  disc,  0  S|t|S  l,and  to  have  zeros 
of  order  ra,  at  t  =  ±  1.   These  integral  representations  provide  the  natur- 
al generalizations  of  the  corresponding  axially  symmetric  representations 
(m  =  0  in  the  above)  that  were  presented  in  an  earlier  work.   As  in  the 
axially- symmetric  case,  they  are  designed  to  satisfy  the  time-reduced 
wave  equation,  the  boundary  conditions,  the  radiation  condition  and  to 


Ill 


have  the  proper  edge  behavior^  assuming  merely  that  the  f's  exist.   A 
straightforward  calculation  reveals  the  f's  to  be  solutions  of  a  Fred- 
holm-type  integral  equation  of  the  second  kind.   Whenever  the  excita- 
tion is  sufficiently  regular  in  a  neighborhood  of  the  aperture^  the 
integral  equations  guarantee,  for  ka  small,  the  existence  of  the  desired 
types  of  f's  and,  moreover,  furnish  approximate  (truncated  Neumann  series)  ex- 
pressions for  them.  Application  of  these  results  is  made  to  the  special 
case  of  plane-wave  excitation.   Approximate  expressions,  in  powers  of  ka, 
are  derived  for  the  far  fields,  transmission  coefficients,  aperture 
fields  and  edge  fields  of  the  solutions  associated  with  the  first  three 
modes  (m  =  0,1,2)  in  the  Fourier-Bessel  expansion  of  the  plane  wave; 
the  leading  term  as  well  as  several  higher  order  terms  are  given  for 
each  quantity.   Analogous  expressions  are  derived  for  m  §  5;  but  here 
only  the  leading  terms  are  given. 
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on  of  the  incident 
e  wave 


A  view  of  the  diffracting  screen.  J     is  the  screen,  (2  is  the  aperture 
of  radius  'a'  centered  at  the  origin  of  coordinates.  J  +   (Z    coincides 
with  the  (x,y) -plane.   The  symbols  v,    <ji,    and  z  are  cylindrical  coordi- 
nates of  an  arbitrary  point  of  space;  r  denotes  the  distance  from  the 
z-axis  (the  axis  of  symmetry);  ^   is  measured  from  the  positive  x-axis 
in  the  (x,y) -plane.   The  angle  0  is  measured  from  the  positive  z-axis 
to  the  ray  that  connects  the  origin  to  (r,^,z).   The  angle  J   fixes  the 
direction  of  propagation  of  the  incident  plane  wave. 


1.   INTRODUCTION 
We  are  concerned  vlth  the  problems  of  diffraction  by  a  circular 
aperture  In  a  perfectly  soft  and  In  a  perfectly  rigid  Infinite,  planar 
screen.   Our  attention  Is  focused  primarily  on  the  lov  frequency  end  of 
the  spectrum  where  the  product  ka  (=a)  of  the  vave  number  k  and  the  aper- 
ture radius  'a'  is  small.   The  solutions  of  J.  Bazer  and  A.  Brown  ^  -'  for 
axl ally- symmetric  excitation  are  herein  generalized  to  cover  the  case  of 
non-axlally  symmetric  excitation  consisting  of  a  superposition  of  'modes' 
of  the  form  cos[m(^-jZ^   )]u    (r,  z),   m  =  0,1, 2,  ...  .   In  this  expression, 
r,  /  and  z  denote  the  cylindrical  coordinates  of  an  arbitrary  space  point, 
with  the  z-axis  serving  as  the  axis  of  symmetry  [see  Figure  l]  ;  j^  ^    is 
an  arbitrary  phase  angle.   As  in  the  axl ally-symmetric  case,  the  new 
solutions  take  the  form  of  Integral  representations,  "^ -.         (f-,    )  and. 
'^P   (-^p    )  ^^y  •>  "^M-Ch  are  related  to  certain  unknown  functions 

f .    (t),  j  =  1, 2  by  meanis  of  single  quadratures.  These  representations 

J 

are  designed  to  satisfy  automatically  the  time-reduced  wave  equation,  the 

boundary  conditions,  the  Sommerfeld  radiation  condition  and  to  have  the 

proper  edge  behavior.    According   as  the  first  or  second  boundary- value 

problem  is  under  consideration,  the  condition  that  the  normal  derivative 

of  the  solution  or  that  the  solution  itself  be  continuous  through  the 

aperture  then  leads  to  Fredholm  integral  equations  of  the  second  kind  for 

for  f^    and  fp    .   Whenever  the  excitation  is  sufficiently  regular  in 
—  .  .  _  _ 

Here,  and  in  what  follows,  it  Is  assumed  that  the  excitation  consists 

of  a  single  mode.   The  superscript  m  identifies  which  harmonic  is  under 
discussion.   Quantities  relating  to  the  soft-screen  problem  or  first  bound- 
arj'-value  problem  are  distinguished  from  their  analogs  in  the  rigid-screen 
or  second  boundary- value  problem  by  appending  the  subscripts  1  and  2. 
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the  neighborhood  of  the  apertiire,  these  integral  equations  guarantee,  for 
sufficiently  small  a,  the  existence  of  f^^'^'s  of  the  desired  type  and, 
moreover,  furnish  approximate  (truncated  Neumann  series)  expressions  for 
these  functions.  We  apply  these  results  to  the  special  case  of  plane- 
wave  excitation.   Specifically,  for  m  =  0,  1  and  2  we  make  use  of  the 
approximate  f^   's  to  obtain  approximate  expressions.  In  powers  of  a,  for 
the  aperture  fields,  far  fields,  edge  fields  and  transmission  coefficients; 
the  leading  terms  as  well  as  higher  order  terms  in  a  are  given.   These  ex- 
pressions evidently  derive  from  the  solutions  associated  with  the  first 
three  hannonics  in  the  Fourier-Bessel  expansion  of  the  plane  wave.  Analo- 
gous expressions  are  derived  for  the  higher  harmonics,  m  g  5^  in  this 
expansion;  but  here  only  the  leading  term  is  given  for  each  of  the  above- 
mentioned  field  characteristics.  When  these  results  are  further  special- 
ized to  the  (axially-symmetric)  case  of  normally-incident  plane  wave  exci- 
tation,  they  are  found  to  agree  with  known  results  obtained  by  other 
means.  Moreover,  it  appears  that  some  of  these  results,  chiefly  those 
associated  with  higher  harmonic  excitation,  are  new. 

The  integral  representations  for  the  higher  harmonics  constitute 
the  essentially  new  element  of  this  paper  from  which  all  other  results 
follow.  As  will  be  seen,  these  representations  can  be  generated  in  a 
natu^ral  way  from  the  axially-symmetric  solutions  presented  in  reference 
[1] .   The  axially-symmetric  solutions  are,  in  turn,  the  analogs  of  simi- 
larly represented  solutions  of  the  first  and  second  boundary- value  prob- 
lems of  potential  theory.  We  refer  the  reader  to  the  book  by  Green  and 


See,  for  example,  C.J.  Bouwkamp's  papers  in  reference  [2j 
See  the  last  paragraphs  of  Sections  5  ^^'^  1- 
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Zerna'-  -J  and  especially  to  the  recent  work  of  W.D.  Collins ^  '  -^    where  simi- 
lar representations  have  been  exploited  to  solve  some  interesting  problems 
in  potential  theory.   Several  authors  have  succeeded  in  reducing  aperture 
and  disc  problems  with  axial  symmetry  and  a  small  to  the  problem  of  solving 

r5i 

Fredholm  integral  equations  of  the  second  kind.   Among  these  are  W.  Magnus  >-  -^ , 

N.  ChakoL^J,  D.S.  Jones '-'''■1 ,  N.I.  Akhiezer  and  A.N.  Akhiezer^-',  B.  Noble  L^-l , 

fiol  r^  cl  * 

R.C.  McCamy  and  A.E.  Heins  ^  -■ ,  and  W.D.  Collins  l  -I  .    Our  integral  equa- 
tions for  the  zeroth  harmonic  are  most  closely  related  to  one  first  obtained 

^j    xj.^.    ^j^i-i^^    .     The  non-axially  symmetric  extensions  of  these  solutions 

*** 
have  not,  to  our  knowledge,  been  published.     Mention  should  here  be  made 

of  J.B.  Keller's  geometrical  theory^  -I  which  requires,  to  be  sure,  that  a 

be  'large',  but  is  limited  neither  to  axi ally- symmetric  excitation  nor  to 

circular  apertures.   It  turns  out  that  the  formulas  for  the  aperture  fields, 

far  fields  and  transmission  coefficients  calculated  in  accordance  with  this 

theory  yield  accurate  results  for  values  of  a  as  low  as  a  =  1.   Most  other 

attacks  on  apertiire  and  disc  problems  have  been  carefully  surveyed  by 

[21 
C .  J .  Bouwkamp  •-  -^  . 

2.   FORMULATION  OF  THE  PROBLEM 
Let  (^  be  a  circular  aperture  of  radius  'a'  in  a  planar  screen  ;a    and 
let  d   +    )/   coincide  with  the  (x,y) -plane,  the  origin  being  located  at  the 
center  of  <^  [see  Figure  l] .   Then  the  solution  U  of  the  first  boundary-value 
problem  may  be  expressed  as  follows  [cf.  reference  [2bJ ,  p.  58] • 


■X- 

This   work  deals   with  the  problem  of  diffraction  by   coaxial  parallel  discs. 

Jones   treats   the   rigid  disc  problem  which,   by  virtue   of  Babinet's  principle, 
is   eqiiivalent  to  the   first  boundary- value  problem  for  the   aperture. 

In  a  private  communication,  W.D.  Collins  has  informed  us  that  he  has  inde- 
pendently arrived  at  our  solution  for  the  case  m  =  1  in  a  study  dealing  with 
the  torsional  oscillations  set  up  in  a  half  space  (elastic)  by  a  disc. 


k  - 


(2.1)     U^  =  U^(r,^z)  = 


U°(r,^z)  -  U°(r,^-z)  +  V^(r,^-z),   z  S  0 


N^{v,4,z), 


z  s  0. 


In  this  equation,  the  triplet  (r,sz(,z)  denotes  the  cylindrical  coordinates 
of  an  arbitrary  point  in  space.   The  angular  coordinate  jZ$  is  measured  from 
the  positive  direction  of  the  x-axis,  as  shown  in  Figure  1,  so  that 


X  =  r  cos 


i^  , 


(2.2) 


y  =  r  s 


In  ^   , 


0  S  0  <  2ir. 


The  incident  wave,  that  is,  U  ,  is  assumed  to  impinge  upon  the  screen  from 
the  left  (in  the  half  space  z  <  O).   The  function  V,  describes  the  scattered 
field;  it  Is  defined  for  all  z  ^  0  and  is  required  to  have  the  following 
properties : 

(i)    V,  is  a  solution  of  the  wave  equation  AV  +  k  V  =  0; 
(11)   V^  =  0  on  /j 

(ill)  V  satisfies  the  Sommerfeld  radiation  condition  at 
infinity} 


(2.5) 


sv,     au"- 


(^^)    ^^  =  ^ 


in^ 


(v)    V,  Is  everywhere  finite; 

(vi)   V  V  is  quadratically  integrable  over  any  domain  in 
3-space,  including  domains  having  some  or  all  points 
of  the  circular  edge  as  boundary  points. 


Similarly,  the  solution  U  (r^jZ^^z)  of  the  second  boundary- value 
problem  may  be  expressed  as  follows. 


r 

U°(r,^z)  +  U°(r,^-z)  -  Vp(r,^-z),    z  S  0 


(2.4)     U2(r,^z)  =  < 


^^{v,i>,7.),  z  g  0. 


Here,  the  scattered  wave  function  Vp  is  defined  in  the  half-space  z  S  0 
and  is  required  to  have  all  the  properties  listed  above  except  (ii)  and 
(iv)   which  must  be  replaced  by 

(ii)  ^^=  0  on     ^  , 

(2.5) 

(iv)'        Vg  =  U°  in  ^    . 

In  both  the  first  and  second  boundary- value  problems,  the  incident 
wave  will  be  assumed  to  have  the  form 


(2.6)     U°(r,^z)  =  cos  m(s^-i!$^^°^^)u^^)(r,z). 


m  =  0,1, 2, . , . 


where  j>   ^^'    denotes  an  arbitraiy  phase  constant  and  \C      (r,  z)  is  a  func- 
tion of  r  and  z  alone.   Moire  general  excitation  can  be  formed  from  a  super- 
position  of  such  'modes'. 


Questions  relating  to  the  convergence  of  this  and  the  corresponding  Fourier 
series  for  the  total  wave  function  U.^  j  =  1,2,  when  an  infinite  number  of 


modes  are  involved^  vlll  be  left  open. 


For  the  purposes  of  the  ensuing  analysis,  it  is  necessai^  to  impose 
suitable  conditions  on  the  incident  wave  at  z  =  0.   To  formulate  these 
conditions,  we  first  write 

(2.7)     u^'"\r,0)  =  r%2^'^^(r)         ,     m  =  0,1,2,..., 


and 

(2.8) 


ou^   (r,  z, 


m   (m)  ,    .  />  T  o 

=  r  u^^  '(r)   ,    m  =  0,1,2^ 

z  io 


for  the  u^^'    and  their  normal  derivatives  in  the  aperture.   In  terms  of 
the  u^^'^'^'s,  k  =  1,2,  m  =  0,1^2,...,  the  aforementioned  conditions  may 
now  be  expressed  as  follows:   It  is  required  (l)  that  the  u^^   (r)'s  be 
even  functions  of  r,  (2)  that  they  be  analytic  functions  of  r  whose 
analytic  continuations  in  the  complex  r-plane  are  regular  at  least  in  a 
(closed)  circular  domain  of  radius  a  +  A  (A  >  0)  centered  at  r  =  0  and 
finally,  (5)  that  they  be  absolutely  integrable  over  every  finite 
r-interval  and  &[e     J   (A  >  O)   for  large  r. 

Consider,  for  example,  the  cases  of  plane-wave  and  spherical-wave 
excitation.   In  cylindrical  coordinates  one  has  the  formulas 

plane  wave: 

/„  ^s  f       /      \  IkTzcosz+rsin/cos^l 

(2.9)  U^(r,0,z)    =   e      L  ^        ^J 

■^^  /j     (krsin7)\ 

ikzcos7      >    .m       ml     m  ^  1      ^      ^ 

m=o  \  r  / 


This  condition  is  imposed  so  that  the  Laplace  transforms  of  the  u^  (r) 
and  of  certain  related  functions  may  exist. 

Time  dependence  of  the  form  exp(-icjt)  is  assumed  throughout,  to  denot- 
ing the  angular  frequency. 


spherical  vave; 


IkR 


(2.19)     U^(r,^z)   =  ^ 


=  2Li  €   cosm(iZ$-s!$    )r 


.    J   (Xr)j   (Xr   )e 
m    /        m  m^      o 


V^2-2 


z-z      i\/X  -k 
o' 


qi=0 


"v^^-^ 


r  */X  -k 


XdX  , 


where    z     <  0, 
o         ' 


(2.11)  € 


m 


and 


1,   m  =  0  , 


2,   m  i.l  ^ 


(2.12)  R  =     yr^+  r  ^-    2rr     cos{<i)-<i>   )   +  (z-z    )^ 


In  the  expressions  for  the  plane  wave,  it  has  been  assumed  that  the  direc- 
tion of  incidence  J  s  in  the  (x,z)-plaiie  and  makes  an  angle  7,0^7  <  Til's., 
with  the   z-axis    [see  Figure   l]  .      In  equation   (2.10)   the  triplet   (r  ,<i)   ,z    ) 

* 

denotes  the   coordinates  oi"  the  point  soiirce. 

Compartng  the   appropriate  hannonics  in  equations    (2.9)    and   (2.10) 
with  the  right  member  of  equation  (2.6)   one  finds,    in  view  of  equations 
(2.7)    and   (2.8),    that 


plane  wave    : 


(m) 


(2.15)        u/°'^(r)    =  ikcosni5^"''(r)   =     i'^^cosre   J   (krsin7)/r'"  = 
1  ^  mm 


*  1dm 

The  last  equation  follows  from  the  recursion  formula  (-  —  — )   Z  fz)  = 
^  ^   z  dz'   o^  ' 

Z  (z)/z  ,  which  holds  for  all  cylindrical  Z   (z)  functions  of  integral 
order  m. 
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ikcos7e   (-1)"^    /,     :,Nm 


(ksinr)"" 


?  £f  Jo(l---in7) 


spherical  vave; 


z_  yJX   -k 


O 


/    \                        p    J   (Xr)j   (Xr   )e 
(2.1M       u,(")(r)  .  -c J  J "      :  ^XdA,  z^<0; 

ij  r 

o 

z    v^  -k 
/    N                       p   J  (Xr)j  (Xr  )e 
(2.15)        u^    '(r)   =     e      /     — — -— —  XdX,  z     <  0    . 

o  r    v^   -k 

The  VL^      (r)'s  in  (2.13)  clearly  have  all  the  properties  listed  above. 
The  same  can  be  shovn  to  hold  true  for  the  u  '  s  of  equations  (2.1^+) 

and  (2.15).   From  this  last  remark  it  follovs  that  our  requirements  on 
the  behavior  of  the  incident  field  entail  no  essential  loss  of  generality; 
for  any  localized  source  may  be  regarded  as  a  superposition  of  spherical 
waves . 

Finally,  it  may  be  mentioned  that  the  problems  of  diffraction  by 
perfectly  soft  and  perfectly  rigid  discs  of  radius  'a'  can  be  solved 
in  terms  of  the  functions  V  and  V^  introduced  in  equations  (2.1)  and 
{2.k) .      In  fact,  it  can  be  shown  that  (see  reference  [2b], p. 59)  that 

U  (r,^z)  -  yJvJ,-z)    ,  z  <  0  , 


(2.16)   U^^^  = 


U^(r,^z)  -  Y^{t,^,z)      ,     z  >  0  , 


and 


'U^(r,^z)   +  V^(r,^-z)    ,  z  <  0 

(2.17)  U^^^  I 

.U^(r,^z)    -  V^(r,^z)      ,  z   >  0. 


where  U,    ,   denotes   the   solution  of  the   soft  or  rigid-disc  problem  accord- 
ing  as   j  =  1  or  2. 
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5-   INTEGRAL  REPRESENTATIONS  OF  THE  SCAITERED  WAVE  V, 

In  reference  [l] ^  where  only  the  case  of  axially-symmetrlc  excita- 
tion was  studied,  the  scattered  wave  V  =  V     was  assiomed  to  have  the  form 

/  ^     p    ikR  (r.zit)    /  ^ 

}        I  ^ikR'^Cr.z^t)    ^ikR"(r,z;t) 
R"^(r,z;t)       R"(r,z}t) 

where  R'*'(r,z;t)    and  R~(r,  z;t)   were   defined  by 


f;L^°^(t)dt  , 


+  /2  2 

(5.2)     R~(r,zjt)  =  y/r     +  (z  ±  iat) 


the  upper  and  lower  signs  on  the  left  corresponding  respectively  to  the 
upper  and  lower  signs  on  the  right.   The  function  f^  '(t)  was  required 
to  be  (l)  an  odd  function  of  t  and  (2),  a  regular  function  of  t  in  the 
circular  domain  defined  by  the -relation  OS  |t|  S  1  +  A/a.   Since  the 
regularity  of  the  incident  wave  U  in  the  domain  Oslr|Sa  +  A^A>0, 
[cf.  requirement  (2)  following  equation  (2.8)  of  the  previous  section] 
ultimately  implied  that  fJ'°^(t)  satisfied  property  (2),  no  generality 
was  lost  in  requiring  that  f^°'^(t)  enjoy  this  property  at  the  very 
outset. 

We  wish  now  to  present  a  generalization  of  the  above  Ansatz.   Owing 
to  our  assumption  about  the  form  of  the  incident  wave  [see   equation  (2.6)J, 


Use  has  already  been  made  of  this  property  in  equation  (5-1) 
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we  may  reduce  the  problem  to  that  of  finding  functions 

(5.3)      V^  =  cos  m((Z^-^^'"))V^^'"V.z)  ,  m=  0,1,2,... 

which  for  each  m  satisfy  conditions  (i),  (ii  ),  (iii),  (v),  and  (vi)  of 
equation  (2.5)  and  the  relation 


av/")(r,z; 


m  (m) /  \ 
=  r  u^  -^(r)  ,      0  s  r  <  a, 

z-l  o 

in  the  aperture  ^.   To  find  such  functions  we  proceed  as  follows.   Let 
V(x,y, z)  be  any  solution  of  the  wave  equation 

(3.5)      4  +  ^-^4-^^'^=0,  z>0. 

dx    by         hz 

Then,  continuing  in  a  purely  formal  manner,  we  find  that  the  functions 
D  V  and  D  V  defined  by 


(5-6)      ^±V  =  (|^  ±  i  |-)v 


are  likewise  solutions.   In  fact,  (D  )'\  defined  by 


I  < m  factors >| 


(5.7)      (D^)\=  i^±^l)    ••••(4*^ 


V 


are  solutions  of  equation  (5-7)  for  each  m,  m  =  1,2,...  .   Moreover,  if 
V  vanishes  on  the  screen,  then  (D_^)  V  will  also  vanish  on  the  screen. 
The  application  of  the  D  's  and  D  's  on  V  may  of  course  produce  unwanted 
singularities  in  the  resulting  wave  function  at  the  circular  edge,  but 
we  shall  ignore  this  possibility  for  the  moment. 
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Representing  D  in  polar  coordinates  [see  equation  (2.2)J,  we  arrive 
at  the  operator  identity 


±±^   ,  b     .    i    h 


(5.8)      D,  =  e--M^ 


*  7^) 


.m 


If  now  (D, )   is  applied  to  any  axl ally- symmetric  solution  -  W  =  W(r, z)  say 
-  of  the  wave  equation,  it  is  found  that 


(5.9)     (Dj\ 


±imjz^ 


(m-1)' 


(m-2) 


W(r,z)  , 


which  may  be  written  more   compactly  as 


/,    ■,^\  /t^    \niT  ±imc^  m/1    a  V 

(5.10)  (D^)^  =  e        '"r  (-^) 


W. 


The  above  considerations  and  this  last  identity  suggest  that  the 
proper  representation  for  the  scattered  field  V  is: 


(5.11)   V^  =  cosm(jZ$-sz^^^'"^)V^^'"^=  cosm(^-jZ^^ 


(m) 


m,l  a  ,m 
^  ^r  ^)  , 


p      ikR"*"(r,z;t)    ,    . 

'      ^  Y"'^(t)dt) 


_^   R"*'(r,z;t) 


In  effect,  we  have  replaced  W  in  equation  (5.9)  ty  an  axially  symmetric 


(o) 


wave  function  having  precisely  the  same  form  as  the  scattered  wave  V, 
of  the  ajclally-symmetric  solution  [cf.  equation  (5.l)j. 

The  functions  fJ''"^(t),  m  =  1,2,...  ,  are  reqiolred  to  share  with 
f^°^(t)  the  property  of  being  (l)  odd  as  functions  of  t  and  (2)  regular 


15  - 


functions  of  t  in  the  clrciolar  domain  defined  by  the  relation  OS  jt]  S  i  +  ^/a 
[see  the  text  following  equation  (5.2)J.   In  addition^  to  ensure  the  proper 
edge  behavior,  it  is  necessary  to  insist  that  for  all  m  S  1 


(5.12) 


dt^ 


=  0  ,       k  =  0,1,  . .  .,ni-l  , 


t=l 


vhlch,  on  account  of  the  oddness  of  f^   (t),  implies  that 


(5-15) 


d\^"^)(t) 

~ S 

dt 


=  0  , 


k  =  0,1, .. . ,m-l 


t=-l 


In  the   above,    as   in  what   follows,    it  is   to  be   understood  that 


;5.1^)  11m     R~(r,z;t)    = 

z;o 


/~2  2~2  ,  2  2.2 

+  Vr 


at      ,    when  r     >  a  t 


This   specification  of  the  branches   of  R   (r, zjt)    immediately  implies   that 


(5.15)  11m     R"(r,z;t) 

Z4.0 


/  2  2  2  ^  ^  . 

+  V  8-  "b     ~   ^     }   when  r     <  a  t 


2  2.2 


At  points  in  the  aperture,  the  only  points  in  physical  space  where 
the  integrand  of  (5.11)  may  be  singular,  V,    and  its  first  and  higher 
order  derivatives  are  defined  to  be  the  limits  approached  when  z  approaches 
zero  through  positive  values.   The  limits  themselves  can  be  derived  with 
the  aid  of  the  properties  of  f    (t)  listed  above  and  Cauchy's  Theorem; 


ih 


Appendix  I  of  reference  [l]  relating  to  the  calculation  of  the  edge  behavior, 
exemplifies  the  sort  of  argtunent  required  .   For  the  sake  of  "brevity,  ve 
shall  not  justify  the  ensuing  limiting  processes.   It  is  enough  to  say 
that  the  several  limiting  processes  employed  in  what  follows  are  justifi- 
able and  agree  with  those  which  are  arrived  at  on  the  basis  of  purely  formal 
considerations . 

Since,  in  Sommerfeld's  notation,  we  have  the  relations 


(5.16)     ^,^'^(v)=7^  H^(^)(w)  = 


iw 
e 

iw 


(5-lT)     ^J'^H=V^  H^^^(^)(w)  , 


(5.18)     (-if)  i   (l)(w)  =  "^      , 
^  ^   w  dw   ^o   ^         m    ' 

w 


and  since  for  each   natural  number  m  we  have  the  identities 


,  2m 
k 


kR"*"  d(kR'*') 


m 


^,^'W"). 


we  may  rewrite  equation   (j-H)    as 


The   requirement  in  property   (2)   that   f^      (t)   be  regular  in  the  domain 
OS    |t|    g  1  +  A/a  is   for  example,    essential  here. 

■x-x-  F"!  p~| 

See  equations  (25)-(55),  pp  ko^-h-o6   in  J. A.  Stratton's  bookL  -■  .   Note 

that  Stratton's  h  ^  "^(w)  corresponds  to  Sommerfeld's  t,   ^      (w)  . 
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(3.20)   V^  =  cosin(sz$-s!^^^"V^^°') 


cosm(  jz^-jz^  ^   ) 


"■      t  ^^)[kR+(r,z;t)] 


^1  [kE-^(r.z;t)]"^    ^ 


-1 


whenever  the  relations  z  >  0,  or  z  g  0  and  r  >  a  h£)ld. 

If  it  Is  assvuned  merely  that  the  f^   's  of  the  desired  type  exists 
then  It  can  be  shown^  even  at  this  stage,  that  the  corresponding  functions, 
V  =  cosin(sz$-jii^  ^™  )V-.    ,  m  =  0,1,...,  satisfy  all  the  conditions  of  equa- 
tion (2.5)  with  the  exception  of  the  aperture  condition  [see  equation  (5«^)] 
That  V-,  satisfies  the  wave  equation  (5-7)  in  the  half-space  z  >0  follows 
directly  from  the  fact  that  (B^)     < exp(ikR  )/r   >  ,  and  linear  combinations 
of  these  functions  enjoy  this  property  for  each  m  g  0.   To  verify  that  V 
vanishes  on  the  screen  [condition  (ii)^  we  need  only  note  that  the  inte- 
grand of  equation  (5.II)  is  non-singular  when  r  >  a  (z  g  O)  and  that  the  • 
f^   's  are  odd  functions  of  tj  the  vanishing  of  the  V-,  's  then  follows 
immediately  on  letting  z  approach  zero.   Next,  we  write 


r  =  R  sin  6  , 

(3.21)      z  =  R  cos  e  ,       0  S  9  ^  rt/2  , 

T,    2    2 

R  =  r  +  z  , 


where   0  is   the   angle  which  the  vector   (x,y, z)    makes  with  the   positive 
z-axls    [see  Figure   ij    and  note   that 

(3.22)  R"(r,z;t)    =   R  ±  iat   cos   0  +  ©tR"""") 
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when  R  is  large.  The  Sommerfeld  radiation  condition  [condition  (iii)^  is 
nov  seen  to  follov  from  this  relation  and  the  known  asymptotic  properties 
of  the  ^  ^  -functions  [see  equations  (3. 19)  a^d.  (5-20)^  or  directly  from 
the  asymptotic  form  of  the  last  member  of  equation  (5-11). 

Since  expTikR  ]/r  and  its  r  and  z-derivatives  are  continuous  func- 
tions of  t  in  the  half-space  z   >  0,   we  conclude  from  equation  (J-H)  that 
for  each  m  g  0  V-  and  grad  V^  are  finite  in  that  region  -  the  singular- 
ity is  not  in  physical  space.   To  determine  the  behavior  of  these  functions 
near  the  edge,  it  is  evidently  enough  to  examine  the  behavior  of  their 
jzi-independent  parts,  namely,  V^  ^   and  grad  V-.    .  We  begin  by  assuming 
that  0  <  |z/a]  <  e  where  e  <  r  is  a  small  positive  number.   Then,  we  de- 
form the  contour  near  t  =  0  in  the  right  member  of  equation  (5-11)  into  a 
semi-circle  of  radius  g  centered  at  t  =  0  and  falling  below  the  real  t- 


ajcis;  lz/a.|  will  then  be  less  than  |t|  everywhere  on  the  contour.   Next, 
we  observe,  that  the  Integral 
follows   [cf.  equation  (j-l)] 


we  observe,  that  the  Integral  representation  of  V^    may  be  expressed  as 


(5.21)  vW  =  Ai|)    <^A  = 


-1 


e 


R'*'(r,z;t) 


ia(z+iat)    ' 


since 


^5-22)     -^  l^^-rr  ia(z+iat)  It  \~r) 

Integrating  by  parts,  we  find,  in  view  of  'edge'  behavior  of  f^^   (t] 
[see  equations  (5-12)  and  (5.15)]  that 


A  circle  has  been   affixed  to  the  integral  sign  to  distingirLsh  between 
the  present  path  of  integration  and  that  used  earlier. 
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(5.23)  v^ 


(m) 


m  .  -.  m-1   p    ikB  (r,  z:t)   ,    ,      ,    . 
_1    ^  (i-^z^t) 


Repeating  this  argiiment  m-1  tlmes^  ve  find  that 


(5.21.)   v^^"^)  =  (^) 


m  n      ^lkR'^(r,z;t) 


^    R"^(r,2:;t)  ^ 


-  (-^) 
dt  ^z+laf* 


-im 


f^^"^)(t)dt 


1 


m 


^  (-4— )!'"  f/°^)(t)dt, 

dt  ^z+iat      1   ^  ' 


where  a  =  ka. 


Now  let  T  be  a  torus  whose  ajcls  of  symmetry  Is  the  z-axis  and 
whose  equation  in  any  cross  section  through  the  z-axis  is  [see  Figure  2j 


(5.26) 


^P  = 


z  =  D  sin  6  , 


r  =  a  +  p  cos  6 


Here,  p  is  the  generating  circle  of  T  and  6  Is  the  angle  measured  clock- 
vise  from  the  screen  to  the  ray  in  the  cross  section  passing  through  the 
edge  and  the  point  P(r, z).  In  these  coordinates  R  becomes 


The  integral  appearing  in  the  first  equation  of  (3. 2k)   has  essentially 
the  same  form  as  that  encountered  in  reference  [ij  at  a  similar  stage  In 
the  development  and  could,  in  fact,  be  manipulated  In  a  similar  manner 
to  obtain  the  edge  behavior.   Here,  we  shall  pursue  a  somewhat  different 
course- 
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i 


,(/t>.8)=P 


■^>Z 


Figure  2 


A  cross -sectional  \rLew  of  the  diffracting  screen. 
The  angle  6  is  measured  from  the  screen  clockwise 
to  the  ray  joining  the  edge  of  the  screen  to  the 
point  P;  p  is  the  distance  from  the  edge  to  the 
point  P. 
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(3.27)    R"*"  =  aVl  -  t^  +  2(p/a)t  e^^  +  (p/a)^  . 


Ignoring  all  but  the  lovest  order  terms  in  (p/a)_,  we  find  that 

,    ,  2    V2     d'^f/'^^(l)  i 

0.28)  V/"')   =     -^  ^     ^  (f)      Bin  6/2  +  a[(p/a)]    . 


ia  dt 


av^^"^         ^    Sv/"^^ 


By  similar  means^    one   can  ascertain  the  behavior  of  — vr —  and  —    - — ^ 
It  turns   out  that 

lav/'")        ^/2    dV^'^V)  i 

(3.29)  (pA)^  -4—  =  -^r2  m         ^W5/2)   +  o[(p/a)2]    , 

ia  dt 


1    ,    av   ("^)  V2     d"^f   ("^^(1)  , 

(3.30)  (p/a)2  i  -1^  =  -^  L_ eos(6/2)    +  ^[(p/a)^]    , 

ia  dt 


which  also  follow  formally  from  equation  (3. 28)  on  applying  the  operators 
S/3p  and  p  h/bb. 

Equations  (3-28)  -  (3'30)  and  the  remarks  made  above  show  that  V, 
and  grad  V,  satisfy  conditions  (v)  and  (vi)  of  equation  (2.3). 

Letting  z  approach  zero  in  equation  (3*11)^  one  finds  with  the  aid 
of  equations  (3-1^)  and  (3-15)  and  the  oddness  of  f^™^(t)  that 


■X- 

One  might  expect  the  last  term  in  equation  (3.25)  to  contribute  a  con- 
stant term  to  this  expression.   However^  it  is  easy  to  show  that  this 
constant  is  zero  owing  to  the  oddness  of  f^   (t).   In  the  second  bound- 
ary value  problem^  however,  the  corresponding  constant  term  does  not 
vanish  [see   equation  (6.3)J. 
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(3.51)    v/"^)  =|£!(i|-)" 

^  ^        1  la       r  dr 


cosh 


[a/vt  - 


(r/a)2j 


(r/a)  Vt^-(r/a)' 


r^(-)(t) 


dt 


In  the  aperture  0  ^   r  <  a.      Moreover^    recalling  that  R  =  R  +  latcosG  +  ©"(R      ) 
one  infers   from  equation   (j.ll)   that 


(5.52)        V^ 


(m) 


.2/      slnh(atcose)f^^"^^(t)dt)     A^  ^)     V" 


r    ^  ikR 


R     -►     00 


where 


(5.55)  A^^°'^(e)   =  -2(iksin9)"'      T    sinh(atcos9)f^^°'^(t)dt 


is  the   far  field  amplitude. 


(m) 


The  transmission  coefficient  t^        ,    associated  with  the  m-th  mode^ 


may  be  expressed  as   follows : 

«/2 


(5.5^)         \ 


(m) 


2rt 

e 

m 


A,^")(G) 


sin  9  d9 


1,      m  =   0 


^        ^m=     < 


2,      m  g   1 

The  numerator  of  this  expression  is  the  energy  flux  in  the  m-th  mode 
that  is  transmitted  to  infinity  on  the  shadow  side  of  the  screen.   The 
denominator  is  energy  flux  in  the  incident  wave  alone  that  is  transmit- 
ted through  the  aperture.   Here,  if  n  denotes  the  outward-directed  normal 
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at  each  point  of  L,    then  the  fliix  of  energy  through  Z  is  taken  to  he 

where  W*  is  the  complex  conjugate  of  W.   From  the  orthogonality  of  the 
set  -Icos   m(jZ^-jZ^    )f  ^  -   0,1,2_,...  ,   it  follovs  that  the  transmission 
coefficient  of  a  sum  of  modes  is  the  sura  of  the  transmission  coefficients 
of  these  modes . 
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k.      IMTEGML  EQUATIONS  FOR  THE  f^^^'s 


There  remains  one  condition  to  be  satisfied  -  the  aperture  condi- 
tion as  formulated  in  equation  (^.k) .      It  will  now  be  shown  that  this 
condition  leads  to  Fredholm  integral  equations  of  the  second  kind  for 
the  unknown  functions  f  J'  '^^  m  =  0^1^2,  .  .  .  .   We  begin  by  calculating 
the  limiting  value  of  hV^^' /bz.      Assuming  that  z  >  0^  we  differentiate 
both  sides  of  equation  (5.11)^  make  use  of  the  identity 


aM[R-(r,z;t)]  .     ^ 

— ^^ =  (z  ±  iat)  -  ^  MrR-(r,z;t)] 


and  the   oddness   of  f  J"      (t)         [m(x)    is   any  differentiable   function  of 
xj ,    and  arrive   at   the   following  result: 


SV^™^  ,      -.    m+1    p       (z+iat)exp[ikR   (r;,z)t)J 


^), 


•J  R    (r,z;t) 


-.    vM-l    p      (z-iat)exp[ikR~(r,z;t)]         ,    , 

r"(i  |:)  /       : f^^")(t)dt. 

^^^  J^  R   (r,z,t)  ^ 


If  we  now  let   Z' -approach  zero,    invoke   the   aperture   condition    [equation 
(3-^)j^    and  make  use   of  the   limiting  behavior  of  R~(r,  z;t)    [see   equa- 
tions   (5.1^)    and   (5.15)3^    ■"■£   find,    in  the   aperture,    0  S   r  <  a,    that 
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r 


{k.2] 


(m) ,    ,         . ,1    h  . 


m+1 


r/a 


exp 


[■VCi"/'  '    ■ 


IcL/Jir/a)    -t 


/y/(r/a)2-t2 


V 


■N 


exp 


.vW  •'  ■" 


ia\/(r/a)    -t 


A/(r/a)2-t2 


tf^^°^)(t)dt 


r 


,1  a    . 
^r  or' 


m+l 


exp 


-o-A/t   -(r/a) 


/a  ^/t^-Cr/a)' 


r/a 


A 


exp 


[aV^ 


(r/a)' 


\/t^-(r/a)2 


tf^^"^)(t)dt 


-,      r/a 
m+l    n         cos 


2i(-^) 
^r  or 


av(r/a)    -t  _ 


/v/(r/a)2-t^ 


tf^^^^(t)dt 


^r  or 


m+l 


■/a  .  r  /; 

sin[aV( 


r/a)2-t2 


2(i^) 
r  dr 


m+l    p        slnh 


A/(r/a)2-t2 
_aVt   -(r/a)  J 


r^^'"^(t 


tf   ^    '(t)dt 


/a  yt2-(r/ 


tf^^")(t)dt. 


r/a 


a)' 


Since   sin 


a.^/{r/a.f-t^\    /  /{ 


O  O  O  Q 

'r/a)    -t   ,    0   S   t     S    (r/a)      is   the   continuation  of 


sinh 


_aVt2-(r/a)2j    /  v/t^-(r/a)2  ,    t^  g    (r/a)^ 


we   conclude   that  for  0  S   r  <  a 


-  ■^'b  - 


(^•3)     (i 


T      r/a 

0 

COS 

a    '^{v/a)    -t  _ 

-  tf^^'")(t)dt 

V(r/a)^-t2 

u^^'^^r)        1  1    ^    °^+l    p      sinh[a  sjt^-{v/a.f_ 

o  V"t   -(r/a) 


tf^^"^)(t)dt 


(o) 


(m), 


Now  define  g^    "^(r)    to  be  u^    '{t)/21l  when  m  =  0^    and 


m-1 


m-2  1 


{k.k)      E^^'^hv) 


21   I    Vl^Vl 


"m-2'^"m-2 


su^    '(s)ds 


m-1 


^     2    a   r^^ 


|a=o 


2"'i(m-l): 


,    2     2.^"-^ 
s(r  -s    )        u 


m-1 


^'")(s)ds  +     2     d 

|J,=  0 


,2^ 


1  m-1 

-^ f      sCl-s^)"""  u.^"^)(rs)ds  +     2     i  ^^^ 

Ifm-l).'       -^  H=o        ^ 


2^(m-l) 


when  m  S  1.   The  d  's  are  constants  to  be  chosen  at  one's  convenience; 

1^ 


*Slnce  u^°'^(s)  Is^  by  hypothesis,  an  even  function  of  s,  it  follows  that 
g  ^"''(r)  is  an  even  function  of  r.   Note  also  that  g^^^^^r)  =  0'[e>cp(Ar)J 
since  u^™^(s)  enjoys  this  property  [see  Section  2}. 
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they  will  be  used  in  Section  5  to  simplify  the  expressions  for  the  g J"™  (r)'s 

associated  with  plane-wave  excitation.   For  the  theory  of  the  present  section, 

one  may  take  them  to  be  zero  or  equlvalently  one  may  absorb  them  into  the 

constants  c  ,g,,...,c   ^  to  be  introduced  presently, 
o  1      m-1 

Operating  on  g^"^' {v)   with  (-  ■^)  >    we  evidently  have 
m 


/I   N      ,1    h  \  (m)  /  s    1    (m)  f    s  ^  T 

(^•5)     (-^)   g^^  '(r)  =  ^  u^^  ^(r),     ms  1 


From  this  equation  and  equation  (^.5)>  it  follows  directly  that 


r/a 
1    d 


'         cos  _a  /\/(r/a)    -t 


(m), 


{k.G)  i^  ,  tf/"'^(t)dt 

r  dr    J  j — - — 2 — 2  "^ 

o  V^r/a)    -t 


./"^(^)  ^  s  ^ 


m-1  „ 

c  2v 

r 


2v+2 
v=:0  a 


■V        p      slnh[a  \/t   -(r/a) 


(m) 


-IJ5?/      :.     . —  tv"'(t)dt 


v/t^-(r/. 


a)2 


The  quantities  c  ,c^...,c    are  at  this  point  arbitrary  constants;  they 
will  be  fixed  later  by  imposing  the  'edge'  conditions  as  formulated  in 
equations  (5.12)  and  (3.15) • 

Equation  (^.6)  is  of  the  same  form  as  the  corresponding  equation 
for  the  axlally-symmetric  case  studied  in  reference  [l]   [cf.  equation 
(^.7)  in  [1]  3  ^^^   ™^y  ^^  treated  in  a  completely  analogous  fashion. 
First,  we  integrate  both  side  of  equation  (^.6)  from  0  to  r  and  obtain 
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the  Integral  equation 


r/a 


cos 


i^.l) 


a 


\f(7/ 


a)    -t 


V(r/a)2-t2 


tf^^'")(t)dt 


m-1 


Sgn 


w 


(s)ds  +     ^^ 


v=o     2(v+l)I 


c  2v+2 


1 


p      slnh|a  -yt   -(r/a) 


Vt^-(r/ 


tf^^"^)(t)dt  -  c^^'") 


a)' 


where 


(k.S)        C. 


(m) 


1 
1 


sinh  as        „  (ra),    ^  , 

sf^        (s)ds 

s  1        ^ 


Next,    we   Introduce   the   variables 


(i^.lO) 


l=t^ 


n  =    (r/a)    , 


into  both  sides  of  equation  {^.J)   which  then  assumes  the  form 


(1^.11) 


r"    ^^i^.W(„ai  =  HW(.)..W(„, 


o     VivT 
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where 


(^.12)  F^^'^ho    =   f/"^(     ^   )/2, 


(^.13)  H^^"^)(t))    = 


^/^\ 


(m) 


m-1  v+1 


1     V 


sg/    ns)ds   +  ^       ^ 


c,^ 


v=o        (v+1) 


(^.IM  j/"^(Tl)=i 


sinh 


[a   VI^ 


vT^ 


F^^'^)(|)d|-C^^"^) 


Employing  the   fact  that  the  Laplace  transform  of  a  convolution  of  two 
functions   is   the  product  of  the  Laplace   transforms^    we   can  solve  equa- 
tion  (^.11)    for  F^        and  express  the   result  as   follows: 


(^.15)  F. 


(m) 


1   d_ 
rt   dT] 


cosh 


[a  V^J 


\Av 


[h/-)(,)+j^^-)(,) 


dn 


,     (m),    ^  .    (m),    V 

=  \        (ti)    +  Ji        (tj) 


Here^    h^    '  {t[)    and  j^      {t\)    are   defined  by 


ik.i6)       h.(-)(n)  =i  ^ 
J-  It 


cosh 


[a  VrPJlJ 


dTl 


A/r\^ 


Hi^°^)(.) 


1  _d_ 
n  dT] 


p-   coshiavTI"!^ 


[a^    ^    ^'^ 


m-1 


a/tR 


g/°'(s)as.2p7;;iTTr^."*')a. 


•i-J  2(v+l) 
v=o       ^  ' 
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/    s  T      ,  p}     cosh  a,   ^Ti-fi  /    N 

o  Vrj-^ 


^.  /^x<"'(^) 


r 


cosh[a   v'Tpii_  sinh  _a.   vT-i^ 


dTl 


•dn 


is/tv^    ^/T■ 


d| 


(m) 


cosh 


(^a   Vt)) 


The   expression  for  h  (tj)   may  be   simplified  somewhat  by  (l)    observing 

that   (-2/a)d[sinh(a  y^t]-h)J    =  [cosh(a   \/ti-^)/    \^t]-[i  Jd|i,    (2)    integrating  by 
parts   and  carrying  out   the  indicated  differentiation^    and  finally  (3) 
introducing  a  new  variabj.e   of  integration,    P,   where   ^i  =  ri  sin  p  .      The 
result  is 


(It. 18)      h/">(n; 


^'^,    ,^" 


rt 


(  m) 
cosh(a   /v^  cosp)sinp  g^    '^(a    /s/t]  sinp)dp 


m-1 


v=o 


«/2 


V      /        V+2 


cosh(a   \fy\  cosp)(sinp)  dp 


With  the   aid  of   Poisson's  formiila  -  i.e. 


2(|) 


n/2 


(i+.l9)        J,(w)  = 


cos(w  cosp)sin     P   dp   ,         t    >  -   i 


v^  r  ( t+i)         o 


-   28 


ve  may  write 


2  r-  Jt/2 

(i^-.20)        h^'^^Cri)   =  /        cosh(a   Vt]  cosp)slnp  g^^°'^(a    V^  sln3)dp 


\'^l  I  ^\ /  r^c^^Ur.       a/^     r.o«RVc;-InR     ;^      ^"^^ 


o 


m-1  J      1  (  iaVn) 


1  Y.     "      -'"^ 


1        V+-2 


2   ^       ^=°        '  (ia   V^  /2)^^2 


Cm)  * 

The   expression  for  j  (t])   may  likewise  be   simplified.      Following  D.S.   Jones   , 

we    (l)    make   use   of  the   addition  formula   for  hyperbolic   functions^    (2)    intro- 
duce obvious    changes   of  variables   of  Integration  and   (3)    carry  out  the   differ- 
entiation with  respect   to  Ti-      The  final  result  is 

,    ,  p         ,,  I  sinhLa(     v^  -    AAi)J        slnhLa(  VT  +    V^) 


-   C   '''"'^    cosh(a   Vn)    /  It 


Combining  equations    (^.15),    (^.20)    and  (^.21),   we   find: 


^2       r-       n/2 
(if. 22)        F^'^^Tl)   =  /        cosh(a   sfy\   cosp)slnp  g^"^^(a    Vn  sinp)dp 


r     1 
See   reference    |_7J    pp   10-11. 
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m-1  J     i(la    V^) 


1         N'  v4      ^-^^ 


+    =  .^        C,,     T] 


2     V^        ^=°        '  (la    v/^/2)^-^* 


1         /■      Cml           )sinh|_a(     vT  -    v^)J        sinh[_a(     //f  +    V^)J  i 
+  —     /   F/"'^(0  < +  — ^>  <i| 


(m) 


C^'"'^    cosh(a    sfv^)  I  n    . 


2  2 

Replacing  t]  by  t     and  ^  by  s     wherever  |   and  y\  appear  in  equation   {\  .22.)    and 

taking  into  account  the  definitions   of  C^^"'      [equation   (^.8)]    and  F^"^   (^) 

[equation   (4-. 12)]    and  the  fact  that   fJ'™'^(t)    is   an   odd  function  of  t,    we   ob- 
tain the   following  Integral  equation  for  the   function  tf^      (t)    for  m  S  1 : 

2  2       "^'"^  im) 

f4   25)        tf  ^™^ft)   =  ^^  "^        [      cosh(atcosp)sinpg^^    '^(atslnp)d3 
1  It        ,/ 


m-1  J      i(lat) 

+     ^^     c     t  — - 

'^        ^=°        "^  (iat/2)^+i 


p       lsinh[a,(t-s)J        cosh(at)sinh(as)  i  ,    . 
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When  m  =  0,  the  siim-term  in  the  right  member  of  this  equation  is  absent. 
Observe  that  the  Bessel  functions  in  equation  (^.25)  are  actually  'ele- 
mentary' transcendental  functions  since 


V+2 


/2  ,  1  d  ^    sin  V 
1     /V  jr  ^  V  dw'     V    -^       }    }    } 


v+2 

w 

From  the  evenness  of  g^    (w)  and  J  i(w)/w  ^,  it  follows  that  the 

1  v+2 

inhomogeneous  terms  of  equation  (^.23)  are  even  functions  of  t.   On  re- 
placing t  by  -t  and  s  by  -s   wherever  s  and  t  appear,  we  find  that 
-tf^  {-t)    satisfies  the  same  integral  equations  as  tf    (t).   Thus, 
assuming  merely  that  the  solution  of  equation  (4.25)  Is  unique,  we  may 
conclude  that  tf^   (t)  is  even  and  hence  that  f^      (t)  is  odd.   Further- 
more,  from  the  behavior  of  g^   (w)   and  J  j^(w)  for  large  w,  we  see 

1  v+2 

that  there  exists  an  N  >  0  such  that  f^  '(t)  =  ©-[e  '  'J  when  |t|  becomes 
large  without  bound  so  that  our  earlier  use  of  the  Laplace  transformation 
is  justified. 

The  above  remarks  are  independent  of  the  particiilar  choice  of  the  c  's 
employed.  We  shall  now  show  how  to  fix  the  c  's  and  how  to  establish  the 
unique  existence  of  the  solution  for  sufficiently  small  values  of  a.   For 
this  purpose,  we  operate  on  both  sides  of  equation  (4.25)  with 
(t   — )  t  ,  n  =  0,1, — ^m-1,   set  t  =  1  and  finally,  enforce  the  'edge' 
conditions  of  equation  (5-12).  We  obtain  as  a  result  the  system  of  equations 


m-1 

(4.24)         2  o^H^^^\a)   =  l^^\ci)   +K^^\a.),        ^i  =  0,1,  . . .  (m-l), 
v=o 


Mr 

See  footnote  on  page  18. 


where 
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(4.25) 


T  ^^)(a) 

V 


H 


2v 


J^_^(iat; 


t=l 


(i^.26) 


I^^^)(a) 


2a      ,1  _d_^^^ 


t/2 


cosh(atcosp)sinPgJ'™^(atsinp)dp 


and 


(^.27)         K^^)(a)  =_i_(iA)^-^ 

/>AF  1    ^  "^^     t 


p       jsinh[a(t-s)]        cosh(at)sinh(ocs) 


t-s 


sfj'")(s)ds 


t=l 


The  elements  of  the  m  X  m  matrix  (t    (an  \i,v  =   0^1^... ^m-1,  are  clearly- 
integral  functions  of  a.   Moreover,,  (^T    (0)J   is  a  triangular  matrix  with 
non-  zero  elements  along  its  main  diagonal  -  in  fact 


(4.29) 


T  ^^^(0)  =  2^vl     V  =  0,1;.  •• -^m-l- 


Thus^  for  sufficiently  small  values  of  |a|^  |a|  <  a     say,T    (a)  has 
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an  Inverse.   This  being  the  case^  we  may,  with  the  aid  of  Cramer's  rule, 
write  the  solution  of  the  system  (^.2^)  In  the  form 


m-1  m-1 

(4.29)     c  =i^^^- ^^ +  ^ 


\l^^\a)  \i^^\cl) 


V 


f 


0   S    \a.\    ^  a,   W^ 
'  '    o 

wherelT  ^'^^(a)|  is  the  determinant  of  (  T  ^^  (a))   and  A  ^*^^(a)  is  the 
cofactor  of  the  element  T  ^'^'^(a).   The  c  's  evidently  Involve  two  types 
of  terms,  those  which  Involve  f^    [through  the  K^^^'s]  and  those  which 

depend  upon  the  gJ-'^'^'s.   If  we  (l)  substitute  the  above  expression  for 

(  m ) 
c  into  equation  (4.23),  (2)  group  the  terms  which  depend  on  f^   with 

the  terms  involving  the  kernel  and  finally,  (3)  group  the  terms  which 

are  independent  of  f^^'   with  the  original  Inhomogeneous  term,  we  arrive 

at  a  new  integral  equation  whose  essential  features  may  be  described  as 

follows:   The  new  Inhomogeneous  term'  shares  with  the  old  the  property  of 

being  a  regular  function  of  t  in  the  disc  OS  [t]  g  1  +  A/a.   In  effect, 

the  domain  of  regularity  of  this  term  is  controlled  by  that  of  g^^^  '^(t) 

[see  footnote  on  page  25^  .   Furthermore,  the  new  kernel  shares  with  the 

original  kernel  the  property  of  being  0(a  )  when  a  is  small.   This  follows 

directly  from  equation  (4.27)  on  representing  the  terms  within  the  curly 

brackets  by  their  power  series.  We  may  therefore  assert  that  the  new 

integral  equation  has  a  unique  continuous  (Neumann  series)  solution  for 
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sufficiently  small  values  of  a.   By  making  a     introduced  above  smaller, 
if  necessary,  ve  may  without  loss  of  generality  assiime  that  the  Neumann 
series  converges  whenever  OS  lal  g  a    .   With  the  unique  existence  of 
a  continuous  f ^ ^   (t)  established  for  0  S  lal  S  a    ,  we  can  infer  from 

1  I   I      o 

equations  ('l.29)  and  (^.25)  -  (k.2"j)    that  the  c  's  are  regular  functions 

of  a,  OS  |a|  S  a  ^'"^  and  that  f  \t)  is  regular  in  the  disc  OS  [t]  S  1  +  A/s 
The  first  remark  follows  directly  from  equations  (^.29)  and  (4.2^)  -  ('+.27). 
To  prove  the  second,  it  is  enough  to  observe  that  (l)  the  new  kernel  is  an 
entire  function  of  t  and  (2)  the  inhomogeneous  term  is  regular  for  all  t  in 
the  disc  OS  |t|  S  1  +  A/a.   Thus  at  least  for  0  S  |a|  S  a  ^'"^  the  unique 
existence  for  the  desired  type  of  f    (t)  is  established. 

In  applications,  it  is  often  the  case  that  the  first  term  in  the  right 
member  of  equation  (4.23)  is  an  analytic  function  of  the  two  complex  vari- 
ables a  and  t  in  the  region  /p  defined  by  0  S  |a|  S  a  ^   , 
OS  |t|  S  1  +  A/a.   Since  the  kernel  is  analytic  for  all  (a, t),  each  term 
in  the  Neumann  series  referred  to  above  is  an  analytic  function  of  (a, t) 
in  ^  .      Moreover,  as  the  convergence  is  uniform  for  (a, t)  in  if  ,    we  may  con- 
clude that  the  limit  function  -  i.e.,  f    (t),  is  analytic  in  the  interior 
of  ^    and  may  be  expressed,  in  this  domain,  as  a  power  series  in  a  with 
power-series  coefficients  in  t.   This  fact  suggests  the  following  scheme 
for  calculating  the  f^       's.   Start  with  equation  (4.23)  and  write  each 
c  as  a  power  series  in  a  -  i.e.,  write 


c  =  2  c  ^^^^"^ 


|i=0 


Next,  set 

00 

fl^"^\t)  =   2  A  ^"^(t)a^ 


u 

|J.=  0 
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►  Z 


Figure    3 


A   cross-sectional  view  of  the   surfaces    (shown   dashed) 

that  boimd  the   domain  oO  .    G        and     (^     are  meridianal 

K       p 

semi-circles  cut  out  of  a  hemisphere  and  torus  whose 
equatorial  planes  coincide  with  the  plane  of  the  screen 
jcl .      Both  torus  and  hemisphere  are  centered  at  0. 
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substitute  this  expression  into  equation  (^.25)  and  equate  like  povers 
of  a-   The  result  is  a  sequence  of  integral  relations  expressing  A    (t) 
in  terms  of  the  c     's  and  integral  terms  involving  A     's  with  lower 
^-indices.   Since  the  kernels  of  the  integral  terms  are  polynomials  in 
s  and  t,  it  follows  that  A    (t)  is  a  polynomial  in  t.   Moreover,  as 
the  first  few  A     's  are  easy  to  calculate  explicitly,  it  follows  that 
the  A    (t)   may  be  calculated  recursively.   The  final  step,  once  the 
A     s  are  known,  is  to  fix  the  c  ^'"^  's  by  imposing  the  'edge'  condi- 

[i  V 

Q 

tions  of  equation  (5.12).   If  terms  of  order  higher  than  a  ,  say,  are 
ignored,  this  procedure  seems  to  be  a  very  efficient  method  for  calcu- 
lating  f    (t).   We  have  made  use  of  it  in  Sections  5  and  7  although, 
in  the  problems,  considered  there,  explicit  formulas  for  the  c  's  are 
easily  found. 

Whenever  a  solution  f^'^'^(t)  of  equation  (^.25)  with  the  required 
properties  exists,  whether  as  a  Neumann  series  or  in  terms  of  the  Fred- 
holm  resolvent,  the  integral  representations  of  Section  5  furnish  a 
solution  V  of  the  aperture  problem  and  equation  (5-11)  gives  the  form 
of  the  edge  behavior.   It  is  easily  shown  that  V,  is  unique  in  the  class 
of  solutions  having  the  properties  listed  earlier  [in  the  neighborhood 
of  equation  (3.3)71  ^^  ^'^'^   radiation  condition  is  strengthened  by  requir- 
ing  k  to  have  a  small  positive  imaginary  part  Im  (k)  >  0.   For,  if  V^'  is 
another  solution  then  v  =  V  -  V  '  enjoys  all  the  properties  of  M.^   except 
that  Sv/5z  vanishes  in  the  aperture.   Consider  the  integral 


/  V(v*-Vv)dW  ; 
here  v^  denotes  the  complex  conjugate  of  v;  <<v  is  the  domain  contained 


For  a  uniqueness  proof  that  does  not  make  use  of  this  added  assumption 
and  which,  in  fact,  weakens  several  others  listed  earlier,  see  the  forth- 
coming report  [13]  by  L.  Levine. 
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within  the  surface  shown  in  Figure  5  (in  cross  section)  and  dW  is  the 
volume  element.   Employing  (l)  Green's  Theorem,  (2)  the  vanishing  of 
V  and  Bv/Sz  on  the  aperture  and  the  screen;,  respectively,  and  (j)  the 
fact  that  V  satisfies  the  wave  equation  and  the  radiation  condition, 
we  find,  on  letting  p  and  R   approach  zero,  that 


/   (|Vv|^-k|v|;dW=0. 


vanishes  everywhere  -  i-e.^  that  V,  -  . 


From  this  and  the  fact  that  Im  (k)  =)=  0,  it  follows  easily  that  v 

This  proof  evidently  applies  when  the  incident  wave  contains  more 
than  one  mode;  however,  the  total  number  of  incident  modes  must  be  fi- 
nite.  When  an  infinite  number  of  modes  are  present,  one  encounters 
problems  relating  to  convergence.   These  require  for  their  resolution 

knowing,  among  other  things,  how  a  '■"'^  and     max    1  f -,    (t)  |  depend 

°        0  S  t  S  1 

on  m  as  m  approaches  infinity.   It  appears  that  complete  success  in 
dealing  with  such  problems  hinges  on  obtaining  simple  explicit  formulas 
for  the  c  's,  possibly  by  enforcing  the  'edge'  conditions  of  equation 
(5.12)  at  an  earlier  stage  and  by  deriving  the  integral  equation  in  a 
different  fashion.   Results  on  this  question  will  be  presented  in  a 
future  report  of  this  series. 

Finally,  it  should  be  mentioned  that  it  is  possible  to  construct 
the  solution  of  equation  (^.25)  in  terms  of  the  eigenfunctions  of  the 
(symmetric)  kernel  sinh  [a(t-s  )] /( t-s)  .   As  D.S.  Jones  L  -I  has  already 
pointed  out,  these  eigenfunctions  are  spheroidal  functions.    However, 


In  a  private  communication  to  J.  Bazer  and  A.  Brown,  CJ.  Bouwkamp  has 
shown  how  to  convert  their  solution  of  the  second  boundary- value  problem 
to  his  original  solution  in  terms  of  spheroidal  wave  functions. 
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it  appears  that  a  solution  in  which  f    (t)  is  so  expressed  would  have 
little  advantage  over  that  obtained  by  employing  spheroidal  wave  func- 
tions at  the  very  beginning. 
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^.   FLAKE  WAVE  EXCITATION  -  THE  FIRST  BOUNDARY  VALUE  PROBLEM 


In  the  present  case,  the  term  involving  g     in  the  integral  equa- 
tion {h.23)    reduces  to  a  simple  expression.   To  derive  this  expression, 
it  seems  that  the  most  efficient  procedure  consists  in  backtracking  some- 
what to  a  suitable  stage  in  the  derivation  of  equation  (^.25).   Let  I     (t)] 
denote  the  first  term  in  the  right  member  of  equation  (^.22).   Then  from 
(i4-.l6),  it  follows  that 


-,'    r»rM:::Vit,^*/T^_iil      /     _.     ""^ 


cosh 


[aA/nniJ 


(5.1)     1/  ^r,)=-  -      sg/  ^(s)ds 


''  ^  r- 

o     y'n-^ 

Choosing  the  arbitrary  constants  which  appear  in  the  definition  of  g    (s) 
[see  equation  (^A)J  in  an  obvious  manner,  we  find,  in  view  of  equation 
(2.13),  that 

I    s  k  cos  y   e  (-1) 

g/  '(s)  =  — -  J  (ks  sm  >). 

^  2(k  sin  r)      ° 

Combining  equations  (5-1)  and  (5.2),  we  then  have 


(m),  ,    '^  "°^  ^  ^m^-"^    d   r   '^"^hCa^ 
(5-2)     V   ^^^  =   ^  ,,   ■   ,,m+l   d^  J    


J-]_(a  Vu  sin  7  )  vTi  in 
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The  -q-dependent  term  in  the  right  member  is  evidently  the  derivative 
of  a  convolution.   Taking  the  Laplace  transform  of  both  sides^  simpli- 
fying the  transforms  of  the  right  member  and  then  taking  inverse  trans- 
forms, we  get 


(5.5)     I^")(n) 


aej-i) 


m 


2it(k  sin  7) 


m 


;inh 


sR 


a  Vn  cos  7 


This,  then,  is  the  first  term  in  the  right  member  of  equation  (^1.22), 

2  2 

which,  on  making  the  substitutions  t^  =  t   and  ^  =  s  ,  becomes 


(5.^; 


tf,^"\t)    = 


/    .  ^m 
ae    (-1) 

t   sinh[at   cos   7] 


jt(k  sin  7  ) 


m 


/—       m-1 


c      t 
2  ^^^  V 

v=o 


2v+2  v+2 

(iat/2) 


jti 


sinhra(t-s)~|    cosh(at)sinh(as )  1      /  s 

_± i )     sf/'"^(s)ds, 

(t-s)  s       ' 


Calculating  tf|°^,  tf^  '    and  tfJ"  '    according  to  the  scheme  outlined 
in  the  preceding  section,  we  find  that 


(5.5) 


(o)  _  (tcos7)a'' 
nk 


t^cos^7   2   2i   3   cos  7    .k  h 
1  +  __ a  +^^     +   -5]—  t  a 


+  1 


2        2 
1    cos  7N     t 

75rt  ^  ^5n  y  +  "T^ 


a5  + 


39 


-h         cos  7    6 
+  —^, —  t 


8lit 


2         71 


a     +Cr(a   )  /       , 


(5.6)  fi 


(1) 


sin27   ,,^4-2^    ^    J     1    ^ 
^  t(l-t    )a       <    5T  + 

rtik 


1            (l+t^Xl+cos^7; 
+  A 

(51)  5: 


a 


+ 


2  k 


i+cosVcos7  n^+^t^^      ^+*    4-  f^^  -  ^±£212.) 
(1+t+t  )  -  yr5T  + ^^2  51       / 


a 


+     Cr(a   ) 


(5-7) 


(2) 


.2cos7sin  7  t(i-t2)2a6  <!^  + 


rtk 


5 


1 120 


(2+cos^7)(t^+2)  1 


5.5; 


a 


+       0~(a-   ) 


Let  e   =   e(r)   be  defined  by 


(5-8) 


V  1  -   (rA 


0  S  r  S   a 


With  the  aid  of  equation  (5-51),  we  obtain  the  following  expressions 
for  V  ^"^^(r.O);  ra  =  0,1^2  in  the  aperture: 


*The  total  aperture  field  is  cosms!^  V^^'{r,0)    (cf.  equation  (5. 11)) 


i+0 


(5.9)      v/°^(r,0)    .I^^S^aa^ 


2         /  2  " 

cos    7      /I        cos    7  \     2 


T-+U  ■  ~T 


2       2i      5 
a     +  -r-  a 
9rt 


a 
■^5: 


k  2  2,^2,2/  8  4  k  2  ' 

cos   7  +  T  cos   7(5-2cos   7)e     +   (1  +  —  cos   7   -  ^  cos   7 


8  2 


^  [-5  "   '°'    V   "■  ^ 


i        2 

e 


a 


6 
cos    7 

7.' 


^    \         f cos    7        2cos    7  \     2 


8lrt' 


720  71 


_^     8    cos    7 


^  2n      , 

cos   7        cos   J  \     4 


5        71  900 


if,        ^2         8  ^  16  6^6 

+  777  (   1  -    2cos    7  +  ^    cos    7   -  -^  cos    7  1  e 


H 


5    —    /    -    55 


a 


+        Q-(a') 


(5-10)      V. 


(1)        2  sin27    ,r.    2k        a 

=  -? (-a.  e  <1  +  — 

3        rt  a  I  "i 


1        cos   7        A 
32  1 2 


cos   y\    2 


a(a^) 


(5.11)    V. 


(2) 


I  •  .2  23  /    "5, 

4icos7sin  7    /£>    ea     +      0"(a   ) 


hi  - 


From  equation  (3-55),  we  obtain  the  following  expressions  of  the 
independent  parts  of  the  far-field  amplitudes: 


(5.12)  a^^^'Hq)  = 


2acos7cos9 
3^ 


a^ll. 


1_ 
10 


2.      2 

cos  0  +  COS  7 


2   2i   3 

97t 


+ 


2    ? 

^^   lOcos  ecos  7      k 
cos  ©  + (-  +   cos  7 


a^  +       &{<x^] 


(5.13)      A^^^^(e)    = 


sin27sin20     4    J  1 


"4  2  2^ 

-^  -  cos    7  +   cos   9 


a 


1 
+   25 


52|  +  |I  cos^r     +  18   cos^r  +  ^  (52  -   28  COBS' )cos' 


COS  0 


a  +   0(a  ) 


(5.1^)   A^^2)(0)  = 


)       .2        2        . 
^acos7sin  7cos0sin  0  6 

1575Tt  "- 


X    a  +i 


2      \  2 

^        cos  7     A        cos   Q 

5  2         r-        2 


2  4 

a     +  Cr(a   ) 


The  corresponding  transmission  coefficients  are: 


kz 


(5.15)      t^(o)   ^  8COS7  ^i.      K        1 


27n 


3  2 

—  +   cos   7 


a 


175 


9       23  2  ^       k 

7j  +  -r-  cos   7  +   :>CQS    7 


a     +       er(a   )  >        , 


(5.16)      t^ 


(l)        32sin  7COS7 
15(^57r)^ 


a       <1  + 


37-21COS   7 


~WT 


a 


283,767  ^     77  2^       253         ^ 


^         ,  6, 
a     +  er(a   ) 


■5.17)      t/2)    ^128cos7sinS        12^^^ 
(I05)(l575n) 


o      2  -, 

2    cos  y 

15-^-9- 


a  +  Cr(a  , 


It  may  be  mentioned,  that  it  is  merely  tedious,  not  difficiilt  to 
obtain  higher  powers  of  a  in  each  of  the  above  expressions  and  to  calcu- 
late the  corresponding  expressions  for  m  =  3-   The  same  applies  to  the 
corresponding  results  listed  in  Section  7- 

In  Appendix  I,  it  is  shown  that 


/  . ^m     .  m  m+1 
/  s      e  (-1)  cos7sin  7a         „     ^ 

(5.18)   f^   (t)  =  (an+l)l:t ^^^-^    ^   ^ 


1  +  e-(a,) 


for  any  m  g  0.   With  the  aid  of  equations  (3.3I);.  (3-53)  and  (5-5^), 
_  . 

The  quantity  €   is  defined  in  equation  (2.11). 


The  order  estimates  in  the  following  expressions  are  conservative. 


45  - 


we   then   find  that 


(5-19)     V. 


(m) 


„m+l      .ra+1    ,         .    .    m 
■d       e   1        mlcos^sm  7 
m 

(2m+l)Ijt 


/  r  ^   .  /-,  ,r^2  m+l 
(-)   Yl-(-)   a 


1  +&{cl) 


(5.20)   a/'^^Cg)  = 


2^     ml  (m+l)  J  a€  (sinGsm?)  cos0cos7  r- 
^       m^  c 
/^  ,  .  >,  ./^  ^s  ■ CL 


2m+2 


(2m+l):(2m+5);rt 


1  +  ©-(a) 


and 


(5-21)   t^ 


(m) 


„6m+6,  ,  N  3r/'   n  ^  iT  3     .2m 
G  2    (ml)   (m+l)l   cos^sm   7  ,   ,,  r 

7^ 7^ ? ct 


rt  [(2mfl):]"[(2nH-3):J 


.l5 


1  +  ©-(a) 


In  the  special  case  of  normal  incidence  (7=0)  our  formulas  for 
V,    (r^O)  and  t,     are  in  agreement  with  those  given  by  C.J.  Bouwkamp 
[see  reference  [2]  pp.  6h-   and  72]  .   We  have  not  succeeded  in  locating 
the  corresponding  expressions  for  "v  ^  0.   Our  results  for  the  higher 
harmonics  m  ^  o  seem  to  be  new.   Observe  that  several  factors  combine 


to  make  V 


^    ^\,     I A    0)  I  and  |t     |  small  when  m  is  large.   First;, 


there  is  the  numerical  factor  involving  the  factorials  of  m.   Next,  there 
is  the  power  of  a  and,,  finally,  ther<?  is  the  factor  involving  powers  of 
sin  7. 
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6.   THE  SECOND  BOIMDAEY- VALUE  PROBLEM 


We  represent  the  scattered  wave  V   (r,^,z)    in  the  half-space  z  >0 


as  follows : 


(6.1)   V^  =  cosm(s;$-s;$^^'"^)V2^'"^(r,z)  = 


cosm(i;i-j!^  ^   ) 


m 


1  a 
r  or 


^   "  f.^")(t)dt 


_^    R+(r,z;t) 


'2 


m  =  0,1,  2,  ..  .  . 

For  each  m,  f    (t)  is  reqiiired  to  be  (l)  an  even  function  of  t,  (2) 
a  regular  function  of  t  in  the  (closed)  circular  domain  0  =  t  S  1  +  A/a 
A  >  0,  and  (5)  to  satisfy  the  'edge'  conditions 


(6.2) 


d%^")(t) 


dt 


t=l 


d^f,^")(t; 


dt 


=  0  ,    k  =  0, 1, . . .,m-l. 


t=-l 


At  points  in  the  aperture,  the  only  points  in  physical  space  where  the 
integrand  can  be  singular,  Vp  and  its  derivatives  are  defined  as  the 
limits  approached  as  z  approaches  zero  through  positive  values. 

Repeating  the  argument  of  Section  3,  we  find  that  conditions  (i), 
(ii)  and  (iii)  [see  equations  (2.3)  and  (2.5)J  are  satisfied.  Apply- 
ing the  method  presented  in  Section  5,  we  find  the  edge  behavior  to  be: 


^5 


(6.5)    V, 


W       ,,  ( 


.^"^(1] 


V"^(l) 


^ill  „    V  111  J  /  _  \ 

m+1      ,^in 
a        at 


(J)   cos(6/2)  +  &[(p/a)]  , 


iS.h)  (p/a) 


1  Sv„ 

2      ^ 


(m) 


dp      m+2 


dt 


m 


cos(6/2)  +  e- 


(p/a)' 


(6.5; 


nH-2 


dt 


m 


sin(5/2)  +  ©•  (p/a)2 


The  variables  p  and  6  are  those  introduced  earlier  in  equation  (5.26) 
fsee  Figure  2].      These  results  imply  that  conditions  (v)  and  (vi)  of 
equation  (2.5)  are  satisfied. 

From  equation  (6.1)  and  equations  (5.1^),  (5-15)^  (5.22),  we  obtain 
the  following  expressions  for  Sv^/Sz  (in  the  aperture)  and  the  far  field 


(6.6)     -^^^ =2cosm(^-^j"^))(i|:) 


cosh 


aVt^-(i-/ 


a)' 


r/a 


Vt^-(r/a)^ 


^   tf2^"^(t)dt  , 


(6.7; 


.(^-^^^"^))A,^™)(9) 


ikR 


V^  ~  cosm(p-^^    ,„2 


R  = 


vK 


where 


(6.8)     A^"^(9)  =  2(iksine)'"   rcosh(atcosG)f2^™''(t)dt  ,   m=   0,1. 


he  - 


The  transmission  coefficient  for  the  m   mode  is 


«/2 


(6.9)     t. 


(m) 


2rt 

6 

m 


A^^™  (0)|  sinGdG 


,  m  =  0,1, 


where  P  is  the  power  in  the  incident  wave  alone  that  is  transmitted 
through  the  aperture  [cf .  Section  Jj  • 

Enforcing  the  aperture  condition,  we  find  for  each  m  g  0,  the  fol- 
lowing analog  of  equation  (^.6): 


(6.10)   i-^ 


1  a 

r 


r/a 
m   p     Cos 


aV(r/a)^-t2. 


V(r/a)2-t2 


f2^"^(t)dt 


a   (m),  .    1  A  S\ 


m  n   sinh 


a 


\/?^ 


a)' 


Vt'-(rA 


f.^"^t)dt. 


Proceeding  as  in  Section  4,  we  are  led  to  the  integral  equation 


(6.11] 


fg^'^t) 


2  _d_. 
It   dt 


_«/2 


t  CO 


sh[atcosp]g^   (atsin(3)sin3d|3 


m-1 


r>    V — '  o     J      i(iat)         1 

+  —    >,  c   t       ,,1  +  ^- 

./-    -*^-J       V  .      ,      V-2  Itl 

-yrt   v=o  /igt.. 


P      sinh[a(t-s)J 


t-s 


f2^"^(t)dt    , 


h7 


where 

(6.12)  g2^°^(r)  =  -^2 


when  m  =  0,  and 


m-1  r  m-1  1 

(6.13)  g/(r)-  S  a^r^-^  f  Vi^V-l  f   V2*V2  '  '  /  ^^g'^''^)*- 


o  o  o 


a        r   /  2   2. m-1   (m  ,  .^ 
— /   s(r  -s  )   u^  '(s)ds 

2"'(in-l) :  J  ^ 

^  o 


2  (m-1):  ^ 


when  m  g  1.   The  d  's  are  constants  to  be  chosen  at  one's  convenience. 
By  arguments  similar  to  those  employed  in  Section  h,    one  can  establish 
the  unique  existence  of  the  desired  type  of  f^   's  for  sufficiently 
small  values  of  a.   In  addition,  the  uniqueness  proof  for  the  total  wave 
function  and  the  scheme  outlined  at  the  end  of  Section  h-   for  calculating 
the  f  ^   ' s  as  power  series  in  a  applies  with  little  change  to  the  pre- 
sent case. 
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7-   PLANE  WAVE  EXCITATION  -  THE  SECOND  BOUNDARY -VALUE  PROBLEM 


Paralleling  the  argument  of  Section  5r  ve  find  the  following  inte- 
gral equation  for  f    (t). 


(7.1)  fg^'^^t)  = 


a(-i)  € 


in-1 


m 


jt(ksin7) 


m 


cosh(atcos7 


>^^S 


2v   Jy-i^i^^^ 


c   t 
V^     ^^       "  (iat/2)''"2 


p      sinh[a(t-s)] 


+   —r 


^.1  (t-s) 


f^^'"n"t)dt   ,  m  =   0,1,  ••• 


Calculating  the   f^        's  by  the  method  outlined  in  the   end  of  Section  k, 
we   obtain  the   following  expressions    for   f^^      ,    fp  and  f^^       : 


(7.2)      f2^°\t)   =1    <ji  -  ^an- 


h        t^ 

J. 

2 


-^  +  ^-   cos    7 


n 


a 


1  +   5cos    78  1  ^2 

Tt 


a 


+ 


2  2        2  !<- 

■2jt     +  1+8  -   2n   cos   7  2,2        cos   7     h 

~Jl ~2^     +^^^ 


1+ 
a  +  ©- 


a 


(7.3)   f2^^)(t)  =  ^'^ 


iiizlkt!)  ,  I  _  3-cosVt^(i^cos^7)   ^2  _^  ^[;3j| 


The  radii  of  convergence  for  these  power  series  in  a  are  not  known. 
CH.  Yang  I-  -^    has  shown  that  it  is  at  least  1/2  in  the  case  m  =  0. 


k9 


(7A)   f2^^^(t)  = 


f^j^ -^       <1  +  &(a) 


,ra 


The  expressions  for  (Sv^   (r,z))/Sz  in  the  aperture  are  accordingly 


Qcf.  equation  (6.1)  and  (6.6)] 


av. 


(o) 


(7-5) 


"^ 


z=o 


2        I,         2i  1 

1     -    CL  -  TT 

itae      )  n  2 


2        2 


2    s    2 


-It    cos    7        /-,    o        2    , 
5 — -  -    (l-2cos   7)e 


a 


>ti 


72-^jt^->t^cos^7   _      2 


3n 


a 


2  22^^  222^ 

52n    -58^+l6it    cos    7-jt    cos    7        8-3fl    cos   7+2^    cos    7      2 


12jt 


F 


6rt 


2      J+ 
5-12cos  7+8cos  7  h 

55 ' 


(7.6) 


av, 


(1) 


z=o 


h   sin  7  r^\  J,   1 


(2cos  7-^)    /o  ),    2^.  2 
-^ 7 +   (2-4cos  7)e 


a^  +  e-(a^) 


Sv, 


(2) 


(7.7; 


z=o 


2     2 
4sin  7  /T,      2  )  ^    a./  \ 
3;^  (-)  a  ^1  +  a(a) 


where  again  e  =  yi  -  (r/a)' 
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The  corresponding  far  field  ampltiudes  are  [see  equation  (6.8)J 


(7.8)  A,(°'(0)=f  ^.4a. 


2   2   2   2 
■2'4-Ht  cos  7-Ht  cos  0 


6/ 


2  2        2  2        2 

23t    -72+5rt   cos   7+5jt    cos   9 


9n 


o?  +     er(a   ) 


(7.9)      k^^\Q] 


-4asin0sin7     2      , ^ 

—  a       <  1   - 


5n 


"I  /'  P  P     '■ 

1    -    Yjr   (|-+COS     7+COS    G) 


a  +  Cf(a   )  >    , 


(7.10)   A2^2)^^)    ^     ^asin^esin^     h     ^^  ^     ^^^^ 


With  the   aid  of  equation   (6-9)    and    these   results,    the   transmission   coef- 
ficients  are    found  to  be 


(7.11)   t, 


(o)  _       8 


2 
n    COS7 


1  +  ^^-?6+>t^cos^7  ^  ^     ^^5) 
9it 


(,.12)    t(l)    =    52sir^^^ 


27n    cos 7 


A9         cos^A     2  ^         ,    ! 


(7.13)   t, 


(2)  128sin   7  8 


a 


50,375it    COS7 


1  +  O-(a^) 
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Applying  the  method  of  Appendix  I;,  we  find^  for  any  integer  m  ^  0, 


that 


/  N       (-1)  e  sm  7a   (1-t  ) 
(7.1^)  f^^'^t)  =  ^ 2 ^ a- 


1  +  (Sr(a 


m+1, 


From  this  expression  it  follows,  as  above,  that 


Sv, 


(m) 


(7.15)      -^ 


-e  m 

m 


i2        (isin7) 


m 


z=o 


(2m):jtae 


m 
,r^        m 
(-)      a 


1  +   ©-(a) 


(7.16)      A2^"^^(0) 


^        mlftn+l^Ie    a(sin7sin9)         „ 
^  m  2m 

a 


(2m);(2m+2):rt 


1  +  (?(a) 


and 


(7.17)   tg^'"^©)  =  (2m): 


^"'^^m:(m+l)l 
(2m);(2m+2).' 


5     .  2m 

V""  ^  ^m1 

^ ^ 

jt  C0S7 


1  +  »(a) 


In  the  special  case  of  normal  incidence  (7=0)  our  formulas  for 
SV^   /Sz,  Ap   (0)  and  t^     are  in  agreement  with  those  given  in 
reference  ["2]  (see  p.  71)-   We  have  not  succeeded  in  locating  the  corres- 
ponding expressions  for  7  /  0.   Our  results  for  the  higher  harmonics  seem 
to  be  new. 


A  misprint  occurs  in  Bouwkamp's  expression  for  A„(0)  which  here  is  de- 
noted  by  A^^   (G)-   It  may  be  rectified  by  replacing  Ap(9)  by  A2(0)A- 
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APPENDIX  I 


We  shall  show  that 


,  V       €  (-1)  a   cos7sin  7t(l-t 

(I.l)     f,^"^(t)  =  ^ a^-^l  1  +  ^(a) 

(2m+l) In 


for  all  integers  m  g  0  when  [see  equation  (2.15)]] 

(1-2)     u/"^(r)  =  i^'-^-^cosre  J  (krsinr)  /r"^. 
-L  mm        / 

We  begin  with  the  observation  that  for  any  A  >  0  and  for  some  suf- 
ficiently small  a_     >  0,  f^   (t)  may  be  expressed  as  a  power  series 
in  a, 


(1.3)      f3_^"^^(t)  =   ^  A  (t)a^  , 


H=o 


which  converges  uniformly  in  t,  0  S  t  g  1  +  (A/a),  whenever  OS  |cl|  S  a 
[see  pp.  55  and  >!+]  •  -^^  explained  on  page  5^,  the  A  's  are  polynomials 
in  t . 

In  the  following,  it  will  be  necessary  to  have  an  estimate  of  the 
degrees  of  A  (t),  A  (t),...,A  (t).   This  is  most  easily  accomplished  by 
expanding  the  right  member  of  equation  ('^.^)     in  powers  of  a;  it  is  then 
easy  to  verify  that  the  degree  of  A  (t),  |a  =  0,1,2,  ...,m  is  at  most  m  +  1. 
But  we  know  that  f^  '^(t)  is  an  odd,  regular  function  of  t  having  zeros  of 


Observe  that  the  regularity  of  J  (w)  for  all  w  implies  that  the  inhomo- 

m 

geneous  term  is  regular  for  all  a  and  t. 
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order  m  at  t  =  +  1  whenever  0  S  |a|  S  a   '' .   Accordingly,  we  may  write 


(l.k)  t^^'^ht)    =   t{l-t^f       ^     B   {t)a}^    ;         0  S  t  S  1  +  A/a  ^ 


\j.=o 


where  the  B  (t)'s  are  constants  (independent  of  t)  when  [i  =   0, 1, 2,  .  .  .  ,m-l 
and  are  even  polynomials  in  t  when  |_i  g  m. 

Now,  combining  equation  C+.J)  and  (1.2)  and  introducing  the  new  vari- 

2  2 

ables  T^  =  (r/a)   and  |  =  t  ,  we  find  that 


(1-5) 


,m+l 


dTl' 


m+1 


^      cos  [a  y^ 


Vnl" 


.  (m) 

"1   (A^  )d| 


sinh 


aVFri 


VFn 


fi^"^VI)d| 


i^ka"^"*"  cos7e  J   {   <^  V^   sin7 


m  m 


„m+l,  /—   ^m 
2    (  Vn  ) 


From  the  form  of  the  left  member  of  this  equation  and  the  fact  that 

J  (  aVri  sinr)  =  i!r(a™),  it  follows  immediately  that  f/"^'^(t)  =  er(a'^)  when 
m  -L 

a  is  small.   Consequently, 


(1-6) 
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in  equation  (l.h).      Moreover,  comparing  coefficients  of  o-  in  both  mem- 
bers of  this  equation,  we  find  that 
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from  which  it  follows  that 
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and  hence  that 
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Combining  this  result  with  equations  (l.6)  and  (l.k),    we  obtain  the 
desired  result,  namely  (1.5)- 
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David  Sarnoff  Research  Center 
Princeton,  New  Jersey 
Attn:  Miss  Fern  Cloak,  Librarian 
Research  Library 

Radio  Corporation  of  America 
Defense  Electronic  Products 
Building  10,  Floor  7 
Camden  2,  New  Jersey 
Attn:  Mr.  Harold  J.  Schrader 

Staff  Engineer,  Organization 

of  Chief  Technical 

Administrator 

The  Ramo-Wooldrldge  Corporation 
P.O.  Box  U51i53  Airport  Station 
Los  Angeles  li5,  California 
Attn:  Margaret  C,  Whitnah, 
Chief  Lllwarian 

Hoover  Microwave  Co. 
9592  Baltimore  Avenue 
College  Park,  Maryland 

Director,  DSAF  Project  RAND 

Via:  Air  Force  Liaison  Office 

The  Rand  Corporation 

1700  Main  Street 

Santa  Monica,  California 

Rantec  Corporation 

Calabasas,  California 

Attn:   Grace  Keener,  Office  Manager 


Northrop  Aircraft,  Inc.  Raytheon  Manufacturing  Company 

Hawthorne,  California  Missile  Systems  Division 

Attn:  Mr.  E,  A.  Freitas,  Library  Dept31ii5  Bedford,  Mass. 

1001  E.  Broadway  Attn:  Mr.  Irving  Goldstein 


Remington  Rand  tJnlv.  -  Division  of  Sperry 

Rand  Corporation 
1900  West  Allegheny  Avenue 
Philadelphia  29,  Pennsylvania 
Attn:  Mr.  John  F.  McCarthy 

R  and  D  Sales  and  Contracts 

North  American  Aviation,  Inc. 
I221U  Lakewood  Boulevard 
Downey,  California 
Attn:  Engineering  Library  U95-11? 

North  American  Aviation,  Inc. 
Los  Angeles  International  Airport 
Los  Angeles  )j5,  California 
Attn:  Engineering  Technical  File 

Page  Communications  Engineers,  Inc. 
710  Fourteenth  Street,  Northwest 
Washington  5,  D.  C, 
Attn:  Librarian 

Phllco  Corporation  Research  Division 

Branch  Library 

U700  Wlssachickon  Avenue 

Philadelphia  hh.    Pa. 

Attn:  Mrs.  Dorothy  S,  Collins 


Raytheon  Manufacturing  Company 
Wayland  Laboratory,  State  Road 
Wayland,  Mass. 
Attn:  Mr.  Robert  Borts 

Raytheon  Manufacturing  Comoany 

Wayland  Laboratory 

Wayland,  Mass. 

Attn:  Miss  Alice  G,  Anderson, 

Librarian 

Republic  Aviation  Corporation 
Farmlngdale,  Long  Island,  N.  Y. 
Attn:  Engineering  Library 

Thru:  Air  Force  Plant  Representative 
Republic  Aviation  Corp. 
Farmlngdale,  Long  Island,  N.Y, 

Rheem  Manufacturing  Company 
9236  East  Hall  Road 
Downey,  California 
Attn:   J.  C.  Joerger 

Trans-Tech,  Inc. 
P.  0.  Box  3U6 
Frederick,  Maryland 


Ryan  Aeronautical  Company 
Lindbergh  Field 
San  Diego  12,  California 
Attn:  Library  -  unclassified 

Sage  Laboratories 
159  Linden  Street 
Wellesley  81,  Mass. 

Sanders  Associates 
95  Canal  Street 
Nashua,  New  Hampshire 
Attn:  N.  R.  Wild,  Ubrary 

Sandla  Corporation,  Sandia  Base 

P.O.  Box  5'^00,  Albuquerque,  New  Mexico 

Attn:  Classified  Document  Division 

Sperry  Gyroscope  Company 

Great  Neck,  Long  Island,  New  York 

Attn:  Florence  W.  Tumbull,  Engr.  Librarian 

Stanford  Research  Institute 

Menlo  Park,  California 

Attni  Library,  Engineering  Division 

Sylvania  Electric  Products,  Inc. 
100  First  Avemie 
Walthall  5U,  Mass. 

Attn:  Charles  A,  Thornhill,  Report  Librarian 
Waltham  Laboratories  Library 

Systems  Laboratories  Corporation 
lli'^52  Ventura  Boulevard 
Sherman  Oaks,  California 
Attn:  Donald  L,  Margerum 

TRG,  Inc. 

17  Union  Square  West 

New  York  3,  N.  Y. 

Attn:  M.  L.  Henderson,  Librarian 

A.  S.  Thomas,  Inc. 

161  Devonshire  Street 

Boston  10,  Mass. 

Attn:   A.  S.  Thomas,  President 

Bell  Telephone  Laboratories 
Murray  Hill 
New  Jersey 

Chu  Associates 
P.  0.  Box  3B7 
Whltcomb  Avenue 
Littleton,  Mass. 

yj-crowave  Associates,  Inc. 
Burlington,  Mass. 

Raytheon  Manufacturing  Company 
Mlssilp  Division 
Hartwell  Road 
Bedford,  Mass. 

Radio  Corporation  of  America 
Aviation  Systems  Laboratory 
225  Crescent  Street 
Waltham,  Mass. 

Lockheed  Aircraft  Corporation 
Missile  Systems  Division  Research  Library 
Box  50U,  Sunnyvale,  California 
Attn:  Miss  Eva  Lou  Robertson, 
Chief  Librarian 

The  Rand  Corporation 
1700  Main  Street 
Santa  Monica,  California 
Attn:  Dr.  W.  C.  Hoffman 

CcHirmander 

AF  Office  of  Scientific  Research 

Air  Research  and  Development  Command 

lUth  Street  and  Constitution  Avenue 

Washington,  D,  C, 

Attn:  Mr,  Ottlng,  SRY 

Westlnghouse  Electric  Corp. 
Electronics  Division 
Friendship  Tnt»l  Airport  Box  7b6 
Baltimore  3,  Maryland 
Attn:  Engineering  Library 
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VOieeler  Laboratories,    Inc. 
122  Cutter  Mill  Road 
Or«at  Neck,    New  York 
Attn:     Mr.   Harold  A.   Wheeler 

Zenith  Plastics  Co. 
Box  91 

Trardena,   California 
Attn:     Mr.   S,   S.   Oleeaky 

Library  Geophysical   Institute 
of  the  TJniversity  of  Alaska 
College 
Alaska 

University  of  California 

Berkeley  b,   California 

Attni     Dr.  Samuel  Silver, 

Prof,   Engineering  Science 
Division  of  Elec.  Eng.   Electronics 
Research  Lab. 

nniversltv  of  California 
Electronics  Research  Lab. 
332  Corv  Hall 
Berkeley  h,   California 
Attn:     J.   R.  Whinnery 

California  Institute  of  Technology 
Jet  Propulsion  Laboratory 
liROO  Oak  Grove  Drive 
Pasadena,    California 
Attn:     Mr.    I.   E.   Newlan 

California  Institute  of  Technology 
1201  E.   Califomi3  Street 
Pasadena,   California 
Attn:      Dr.   C.    Papas 

Carnegie  Institute  of  Technology, 
Schenlev  Park 

Pittsburgh  13,    Pennsylvania 
Attn:      Prof.  A.   E.   Heins 

Cornell  University 

School  of  Electrical  Engineering 

Ithaca,    New  York 

Attn:     Prof.  G.   C,   Dalman 

I'niversity  of  Florida 

Department  of  Electrical  Engineering 

Gainesville,   Florida 

Attn:      Prof.   M.   H.   Latour,    Library 

Library 

Georgia  Institute  of  Technology 

Engineering  Experiment  Station 

Atlanta,  Georgia 

Attn:  Mrs.  J.H.  Crosland,  Librarian 

Harvard  University 

Technical  Reports  Collection 

Gordon  McK'ay  Library,  303A  Pierce  Hall 

Oxford  Street,  Cambridge  3f^,  Mass. 

Attn:  Mrs.  E.L.  Hufschmidt,  Librarian 

Harvard  College  Observatory 
60  Garden  Street 
Cambridge  3^,  Mass. 
Attn:  Dr.  Fred  L.  Whipple 

Hniversity  of  Illinois 
Documenta  Division  Library 
Urbana,  niinois 

University  of  Illinois 
College  of  Engineering 
Urbana,  Illinois 

Attn:   Dr.  P.  E.  Moyes,  Department  of 
Electrical  Engineering 

The  Johns  Hopkins  University 
Homewood  Campus 
Department  of  Physics 
Baltimore  1*5,  Maryland 
Attn:  Dr.  Donald  E.  Kerr 

Sandia  Corporation 
Attn:  ^Organization  l'i23 
Sandia  Base 
Albu-njernue,  New  Mexico 


Applied  Physics  Laboratory 
The  Joh^ Hopkins  University 
B621  Georgia  Avenue 
Silver  Spring,  Maryland 
Attn:  Mr.  George  L.  Seielstad 

Massachusetts  Institute  of  Technologv 

Research  Laboratory  of  Electronics 

Room  20B-221 

Caniiridge  39,    Massachusetts 

Attn:     John  H.   Hewitt 

Massachusetts   Institute  of  Technology 

Lincoln  Laboratory 

P.  0.   Box  73 

Lexington  73,  Mass, 

Htn:   Document  Roan  A-229 

University  of  Michigan 
Electronic  Defense  Group 
Engineering  Research  Institute 
Ann  Arbor,  Michigan 
Attn:   J.  A,  Boyd,  Supervisor 

University  of  Michigan 
Engineering  Research  Institute 
Radiation  Laboratory 
Attn:   Prof.  K.  M.  Siegel 
912  N.  Main  St., 
Ann  Arbor,  Michigan 

University  of  Michigan 
Engineering  Research  Institute 
Willow  Run  Laboratories 
Willow  Run  Airport 
Ypgilanti,  Michigan 
Attn:  Librarian 

University  of  Minnesota 

Minneapolis  lb,  Minnesota 

Attn:  Mr,  Robert  H,  Sturon,  Library 

Northwestern  University 
Microwave  Laboratories 
Evanston,  Illinois 
Attn:   R.  E.  Beam 

Ohio  State  University  Research  Found. 
Ohio  State  University 
Columbus  10,  Ohio 
Attn:  Dr.  T.S.  Tice 

Dept.  of  Elec.  Engineering 

The  University  of  Oklahoma 
Research  Institute 
Norman,  Oklahoma 

Attn:  Prof.  C,  L.  Farrar,  Chairman 
Electrical  Engineering 

Polytechnic  Institute  of  BrooV^lyn 
Microwave  Research  Institute 

?5  Johnson  Street 

Brooklyn,  New  York 

Attn:  Dr.  Arthur  A,  Oliner 

Polytechnic  Institute  of  Brooklyn 
Microwave  Research  Institute 
??  Johnson  Street 
Brooklyn,  New  York 
Attn:  Mr.  A,  E,  Laemmel 

Syracuse  University  Research  Institute 
Collendale  Campus 
Syracuse  10,  New  York 
Attn:  Dr.  C.  S.  Grove,  Jr. 

Director  of  Engineering  Research 

The  University  of  Texas 
Elec.  Engineering  Research  Laboratory 
P.  0.  Box  8026,  University  Station 
Austin  12,  Texas 
Attn:  Mr,  John  R,  Gerhardt 
Assistant  Director 

The  University  of  Texas 

Defense  Research  Laboratory 

Austin,  Texas 

Attn:  Claude  W.  Horton,  Physics  Library 

University  of  Toronto 

Department  of  Electrical  Engineering 

Toronto,  Canada 

Attni   Prof,  G.  Sinclair 


Lowell  Technological  Institute 
Research  Foundation 
P.  0.  Box  709,  Lowell,  Mass, 
Attn:  Dr.  Charles  R.  Mingins 

University  of  Washington 

Department  of  Electrical  Engineering 

Seattle  5,  Washington 

Attn:  G.  Held,  Associate  Professor 

Stanford  University 
Stanford,  California 
Attn:  Dr.  Chodorow 

Microwave  Laboratory 

Physical  Science  Laboratory 

New  Mexico  College  of  Agriculture 

and  Mechanic  Arts 

State  College,  New  Mexico 

Attn:   Mr.  H.  W.  Haas 

Brown  University 
Department  of  Electrical  Engineering 
Providence,  Rhode  Island 
Attn:  Dr.  C.  M.  Angulo 

Case  Institute  of  Technology 

Cleveland,  Ohio 

Attn:   Prof.  S,  Seeley 

Columbia  University 

Department  of  Electrical  Engineering 

Momingside  Heights 

New  York,  N.  Y. 

^ttn:     Dr.  Schlesinger 

McGill  University 

Montreal,  Canada 

Attn:  Prof.  G,  A.  Woonton 

Director,  The  Eaton  Electronics 

Research  Lab. 

Purdue  University 

Department  of  Electrical  Engineering 

Lafayette,  Indiana 

Attn:  Dr.  Schulta 

The  Pennsylvania  State  University 

Department  of  Electrical  Engineering 
University  Park,  Pennsylvania 

University  of  Pennsylvania 

Institute  of  Conperative  Research 

IliOO  Walnut  Street 

Philadelphia,  Pennsylvania 

Attn:  Dept.  of  Electrical  Engineering 

University  of  Tennessee 
Ferris  Hall 
W.  Cumberland  Avenue 
Knoxville  16,  Tennessee 

University  of  Wisconsin 
Department  of  Electrical  Engineering 
Madison,  Wisconsin 
Attn:  Dr.  Scheibe 

University  of  Seattle 

Depjartment  of  Electrical  Engineering 

Seattle,  Washington 

Attn:  Dr.  D.  K.  Reynolds 

Wayne  University 

Detroit,  Michigan 

Attn:   Prof.  A.  F.  Stevenson 

Electronics  Research  Laboratory 
Illinois  Institute  of  Technology 
3300  So,  Federal  Street 
Chicago  16,  Illinois 
Attn:  Dr.  Lester  C,  Peach 
Research  Engineer 

Advisory  Group  on  Electronic  Parts 

Room  103 

Moore  School  Building 

200  South  33rd  Street 

Philadelphia  b,  Pennsylvania 
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Ionosphere  Research  Laboratory 

Pennsylvania  State  College 

State  College,  Pennsylvania 

ATTNi   Professor  A.  H,  Waynick,  Director 


Institute  of  Mathematical  Sciences 

25  Waverly  Place 

New  York  3,  New  York 

ATTNi  Librarian  O) 

Electronics  Diviaion 
Rand  Corporation 
1700  Main  Street 
Santa  Monica,  California 
ATTNi  Dr.  Robert  Kalaba 

National  Biareau  of  Standards 

Washington,  D.  G, 

ATTNi  Dr.  W,  K.  Saunders 

Applied  Mathematics  and  Statistics  Lab. 

Stanford  University 

Stanford,  Qilifomla 

ATTN:  Dr.  Albert  H.  Bowker 

Departmert  of  Physics  and  Astronomy 
Michigan  State  College 
East  Lansing,  Michigan 
ATTN:  Dr.  A.  Leitner 

UnlversitT  of  Tennessee 
Knoxville,  Tennessee 
ATTN:  Dr.  Fred  A.  Flcken 

Lebanon  Valley  College 
Annvllle,  Pennsylvania 
ATTN:   Professor  B.H.  Bissinger 

General  Atomic 

P.  0.  Box  608 

San  Diego  12,   California 

Attn:      Dr.    Edward  Gerjuoy 

Department  of  Physics 
Amherst  College 
Amherst,   Mass. 

ATTN:     Dr.  Arnold  Arons 

California  Institute  of     Technology 
1201  E.  California  Street 
Pasadena,   California 
ATTN:     Dr.   A.  Erdelyl 

Mathematics  Depertment 

Stanford  University 
Stanford,  California 
ATTN:     Dr.   Harold  Levlne 

University  of  Minnesota 

Mlnnpapolls  I'a,   Minnesota 

ATTNi     Professor  Paul  C.   Rosenbloom 

Department  of  Mathematics 

Stanford  University 

Stanford,   California 

ATTN:      Professor  Bernard  Epstein 

Applied  Physics  Laboratory 
The  JohrP Hopkins  University 
8621  Georgia  Avenue 
Silver  Spring.   Maryland 
ATTN:     Dr.  B,   S.   Gourary 


(2)Exchange  and  Gift  Division 
The  Library  of  Congress 
Washington  25,    D.  C, 

Electrical  Engineering  Department 
Massachusetts  Institute  of  Technology 
Cambridge  39,    Mass. 
ATTN:     Dr.  L.   J.  Chu 

Nuclear  Developnent  Associates,    Inc. 

5  New  Street 

White  Plains,    New  York 

ATTN:      Library 

California  Institute  of  Technology 

Electrical  Engineering 

Pasadeni,   California 

ATTN:     Dr.   Zohrab  A.   Kaprielian 


Dr.   Rodman  Doll 
311  W.   Cross  Street 
Tpsllanti,   Michigan 

California  Inst,    of  Technology 
Pasadena,    California 
ATTN:     Mr.  Calvin  Wilcox 

Mr.   Robert  Brockhurst 

Woods  Hole  Oceanographic   Institute 

Woods  Hole,   Mass. 

National  Bureau  of  Stardards 
Boulder,   Colorado 
ATTN:     Dr.   R.  Gallet 

Dr.   Solomon  L.   Schwebel 

3689  Louis  Road 

Palo  Alto,     Salifomia 

University  of  Minnesota 
The  University  of  Library 
Minneapolis  Ih,    Minnesota 
ATTN:     Exchange  Division 

Department  of  Mathematics 
Oniveraitv  of  Califomla 
Berkeley,    California 
ATTN:      Professor  Bernard  Friedman 

Lincoln  Laboratory 

Massachusetts   Institute  of  Technology 

P.   0.  Box  73 

Lexington  73,   Massachusetts 

ATTN:     Dr.   Shou  Chin  Wang,    Room  C-351 

Melpar,  Inc., 

3000  Arlington  Boulevard 

Falls  Church,  Virginia 

ATTN:  Hr.  K.  S.  Kelleher,  Section  Head 

Hq.  Air  Force  Cambridge  Research  Center 

Laurence  G,  Hanscom  Field 

Bedford,  Mass. 

ATTN:  Mr.  Francis  J.  Zucker,  CRRD 

Hq.  Air  Force  Cambridge  Research  Center 

Laurence  G.  Hanscom  Field 

Bedford,  Mass. 

ATTNt  Dr.  Philip  Newman,  CRRK 


Mr.  N,  C.  Gerson 

Trapelo  Road 

South  Lincoln,  Mass, 

Dr.  Richard  B.  Barrar 
Systems  Developmient  Corp. 
2U00  Colorado  Avenue 
Santa  Monica,  California 

Columbia  University  Hudson  Laboratories 
P.O.  Box  239 

lli5  Palisade  Street,  Dobbs  Ferry,  N.  Y. 
ATTN:  Dr.  N.  W.  Johnson 

Institute  of  Fluid  Dynamics 
and  Applied  Mathematics 
University  of  Maryland 
College  Park,  Maryland 
ATTN:  Dr.  ElUott  Monrtroll 

Department  of  Electrical  Engineering 

Washington  University 

Saint  Louis  5,  Mo. 

ATTN:   Professor  J.  Van  Bladel 

Department  of  the  Navy 

Office  of  Naval  Research  Branch  Office 

1030  E.  Green  Street 

Pasadena  1,  California 

Brandels  University 
Waltham,  Mass. 
ATTN:  Ubrary 

General  Electric  Company 
Hlctowave  Laboratory 
Electronics  Division 
Stanford  Industrial  Park 
Palo  Alto,  California 
ATTN :  Library 


Smyth  Research  Associates 
3555  Aero  Court 
San  Diego  3,  California 
ATTN:  Dr.  John  B,  Smyth 

Electrical  Engineering 
California  Institute  of  Technology 
Pasadena,  California 
ATTN:  DrXieorgea  G»  Weill 

Naval  Research  Laboratory 

Washington  25,  D,  C. 

ATTN:  Henry  J,  Passerinl,  Code  5278A 

Dr.  George  Kear 
5  Culver  Court 
Orlnda,  California 

Brooklyn  Polytechnic 
85  Livingston  Street 
Brooklyn,  New  York 

ATTN:  Dr.  Nathan  Marcuvltz 

Department  of  Electrical  Engineering 
Brooklyn  Polytechnic 
85  Livingston  Street 
Brooklyn,  New  York 
ATTN:  Dr.  Jerry  Shmoys 

Department  of  Mathematics 
University  of  New  Mexico 
Albuquerque,  New  Mexico 
ATTN:  Dr.  I.  Kol^dner 

Mathematics  Department 
Polytechnic  Institute  of  Brooklyn 
Johnson  and  Jay  Street 
Brooklyn,  New  York 
ATTNi  Dr.  Rarry  Hochstadt 

Ballistics  Research  Laboratory 
Aberdeen  Proving  Grounds 
Aberdeen,  Maryland 
ATTNi  Dr.  Pullen  Keats 

Dr.  Lester  Kraus 
U935  Whitehaven  Way 
San  Diego,  Califomla 

University  of  Minnesota 
Institute  of  Technology 
Minneapolis,  Minnesota 
Attn:  Dean  Athelston  Spilhaus 

Ohio  State  University 
Columbus,  Ohio 
Attn:  Prof.  C.  T.  Tai 
Department  of  Electrical  Eng. 

Naval  Research  Laboratories 

Washington  25,  D.  C. 

Attn:  W.  S.  Ament,  Code  5271 

Naval  Research  Laboratory 
Washington  2$,    D.  C. 
Attn:  Dr.  Leslie  G.  McCracken,  Jr. 
Code  3933A 

Office  of  Naval  Research 
Department  of  the  Navy 
Attn:  Geophysics  Branch,  Code  bl6 
Washington  25,  D.  C. 

Office  of  Chief  Signal  Officer 
Signal  Plana  and  Operations  Division 
Attn:  SIGOL-2,  Room  20 
Com.  Liaison  Br.,  Radio  Prop.  Sect. 
The  Pentagon,  Washington  25,  D.  C. 

Defence  Research  Member 
Canadian  Joint  Staff 
2001  Connecticut  Street 
Washington,  D.  C. 

Central  Radio  Prop.  Lab, 
National  Bureau  of  Standards 
Attn:  Technloal  Reports  Library 
Boulder,  Colorado 

U.  S.  Weather  Bureau 
U.  S,  Department  of  Conanerce 
Washington  25,  D.  C, 
Attni  Dr.  Harry  Wexler 
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Federal  Comiminlcationa  Commlaslon 

WashltKton  25,    D.  C. 

Attnt     Mrs.   Barbara  C.   Grimes,   Librarian 

Upper  Atmosphere  Research  Section 
Central  Radio  Propagation  laboratory 
National  Bureau  of  Standards 
Boulder,   Colorado 


Argonne  National  Laboratory 

P.O.  Bo3t  299 

Lemont,  Illinois 

Attn:  Dr.  Hoylande  D.  Young 

Bell  Telephone  Labs. 

Murray  Hill,  New  Jersey 

Attni  Dr.  S.  0.  Rico,  3B  -  203 

Carnegie  Institute  of  Washington 
Dept.  of  Terrestrial  Magnetism 
52I1I  Broad  Branch  Road,  N.  W, 
Washington  15,  D.  C. 
Attn:  Library 

Georgia  Tech  Research  Institute 
225  N.  Avenue,  N.  W. 
Attn:  Dr.  James  E.  Boyd 
Atlanta,  Georgia 

nniversity  of  Maryland 
College  Park,  Maryland 
Attn:  Dr.  A.  Weinstein 
Institute  of  Fluid  Dynamics 

Massachusetts  Institute  of  Technology 

Lincoln  Laboratory 

P.  0.  Box  73 

Lexington  73,  Massachusetts 

Attn:  Prof.  Radford,  Division  3  Head 

Willow  Run  Research  Center 
University  of  Michigan 
Willow  Run  Airport 
Tpsllanti,  Michigan 
Attn:  Dr.  C.  L.  Dolph 

School  of  Engineering 
New  York  University 
University  Heights 
New  York,  New  York 

Shell  Fellowship  Committee  of  the 
Shell  Companies  Foundation,  Inc. 
50  West  50th  Street 
New  York  20,  N.  Y. 
Attn:  Mr.  J.  R.  Janssen 

Esao  Research  and  Engineering  Co. 

P.  0.  Box  51 

Linden,  New  Jersey 

Attn;  Mr.  C.  L.  Brown,  Manager 

Union  Carbide  and  Carbon  Corp. 
30  E.  U2nd  Street 
New  York  17,  New  York 
Attn:  Mr.  L.  E.  Erlandson 

Convalr 

San  Diego  12,   California 

Attn:     Mr.  Marvin  Stem 

Bell  Telephone  Labs.,    Inc. 

Ii63  West  Street 

New  York  13,    N.   Y. 

Attn;     Dr.   Mervin  J.   Kelly 

Engineering  Library 
University  of  California 
I1O5  Hllgard  Avenue 
Los  Angeles  2U,   California 

Convalr,   A  Division  of  General  Dynamics  Corp. 

Daingerfield,   Texas 

Attn:     J.   E.   Arnold,    Division  Manager 

Convalr,    A  Division  of  General  Dynamics  Corp. 

San  Diego  12,   California 

Attn:     R.   L.  Bayless,   Chief  Engineer 

Convalr,    a  Division  of  General  Dynamics  Corp. 

San  Diego  12,   California 

Attn:     K.   J.  Bossart,   Chief  Engineer-WS107A 


Convalr,   A  Dlv.   of  General  T^amics 

Corp. 
Fort  Worth  1,   Texas 
Attn:     F.   W.   Davis,   Chief  Engineer 

Convalr,   A  Dlv.   of  General  Dynamics 

Corp. 
Pomona,   California 
Attn:     C.   D.    Perrlne 

Ass 't  Dlv,   Manager,   Engln. 

Shell  Developnent  Company 

Exploration  arri  Production  Res.   Dlv. 

3737  Bellalre  Boulevard 

Houston  25,   Texas 

Attn;     Miss  Aphrodite  Mamoulliles 

RCA  Laboratories 
Princeton,   New  Jersey 
Attn:     Dr.  Charles  Polk 

Stanford  Research   Institute 
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